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Abstract 



In the first chapter we give an introduction to hard diffractive scattering in QCD 
to introduce basic concepts and terminology, thus setting the stage for the following 
chapters. 

In the second chapter we make predictions for nondiagonal parton distributions 
in a proton in the LLA. We calculate the DGLAP-type evolution kernels in the 
LLA, solve the nondiagonal GLAP evolution equations with a modified version of 
the CTEQ-package and comment on the range of applicability of the LLA in the 
asymmetric regime. We show that the nondiagonal gluon distribution g(xi, X2, t, /i 2 ) 
can be well approximated at small x by the conventional gluon density xG(x,fi 2 ). 

In the third chapter, we discuss the algorithms used in the LO evolution program 
for nondiagonal parton distributions in the DGLAP region and discuss the stability 
of the code. Furthermore, we demonstrate that we can reproduce the case of the 
LO diagonal evolution within less than 0.5% of the original code as developed by 
the CTEQ-collaboration. 

In chapter 4, we show that factorization holds for the deeply virtual Compton 
scattering amplitude in QCD, up to power suppressed terms, to all orders in per- 
turbation theory. Furthermore, we show that the virtuality of the produced photon 
does not influence the general theorem. 
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In chapter 5, we demonstrate that perturbative QCD allows one to calculate 
the absolute cross section of diffract ive exclusive production of photons at large 
Q 2 at HERA, while the aligned jet model allows one to estimate the cross section 
for intermediate Q 2 ~ 2GeV 2 . Furthermore, we find that the imaginary part of the 
amplitude for the production of real photons is larger than the imaginary part of the 
corresponding DIS amplitude, leading to predictions of a significant counting rate 
for the current generation of experiments at HERA. We also find a large azimuthal 
angle asymmetry in ep scattering for HERA kinematics which allows one to directly 
measure the real part of the DVCS amplitude and hence the nondiagonal parton 
distributions. 

In the last chapter, we propose a new methodology of gaining shape fits to 
nondiagonal parton distributions and, for the first time, to determine the ratio r\ 
of the real to imaginary part of the DIS amplitude. We do this by using several 
recent fits to F 2 (x,Q 2 ) to compute the asymmetry A for the combined DVCS and 
Bethe-Heitler cross section. The asymmetry A, isolates the interference term of 
DVCS and Bethe-Heitler in the total cross section, in other words, by isolating the 
real part of the DVCS amplitude through this asymmetry one has access to the 
nondiagonal parton distributions for the first time. Comparing the predictions for 
A against experiment would allow one to make a prediction of the shape, though 
not absolute value, of nondiagonal parton distributions. 

In the appendix, to illustrate an application of distributional methods as dis- 
cussed in chapter 4, we will show, with the aid of simple examples, how to make 
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simple estimates of the sizes of higher-order Feynman graphs. Our methods enable 
appropriate values of renormalization and factorization scales to be made. They 
allow the diagnosis of the source of unusually large corrections that are in need of 
resummation. 
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CHAPTER 1 
Introduction 



1.1. What is Diffraction in QCD ? 

The first question we will attempt to answer is: What is diffraction in QCD? The 
answer is, of course, not 42 [[jj but rather complicated and multi layered going to 
the heart of our understanding of QCD or lack thereof. 

Before proceeding, we will introduce customary notation for variables used in 
describing inelastic and diffractive phenomena. 

In the process 

e(l)+p(p) -> e'(l')+X + P'(p') 
or 

P(l)+P(p) P(1')+P(P') 
or 

P(1)+P(P) - X+ V ip') (1.1) 

with the momenta of the particles given in brackets. For defineteness and since this 
type of process will be of greatest relevance in the remainder of the thesis, we will 
introduce the relevant kinematics for the electron proton scattering. The center of 
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mass energy squared is s = (I + p) 2 and the center of mass squared of the hadronic 
system is 

W 2 = (q + p) 2 = -Q 2 + 2p ■ q + m 2 p = Q 2 (-1 + ~) + m 2 p , (1.2) 

with q = (—xp + ,Q 2 /2xp + ,0±) [] being the four momentum of the virtual photon 
exchanged between the lepton and the rest of the system (see Fig. |TTT|) and Q 2 = —q 2 . 
The Bjorken x is defined as: 



o_ 2 

2p ■ q 



x bj = 77— (1-3) 



with p being the four momentum of the incoming proton. 



ep scattering, but also the rest of the reactions in ( |1.1| ), is said to be in the 
diffractive region if X is sufficiently separated from the final state proton; the pre- 
cise meaning of this statement will be explained below. X is used in diffractive 
interactions to denote the part of the final state which might be in the case of ep 
scattering, a real photon, a vector meson or two jets. If the final state proton dis- 
sociates as well, we will denote this state by Y to indicate that it is still separated 
from the state X. The kinematics are the following: The square of the momentum 
transfer at the proton vertex is 

t=(p-p') 2 (1.4) 

where p' is the four momentum of the outgoing proton. In the case of dissociation 
1 We define a vector in light cone coordinates by: 



, v° -I- V 3 V° — V 3 



2 



♦ 
♦ 



♦ 



♦ 




p 



p 



Figure 1.1: Typical lepton-proton scattering with a final state X clearly separated 
from the final state proton with the exchanged object between the virtual photon 
and the proton denoted by a dashed line. 

of the final state proton the proton momentum p' should be replaced by the four 
momentum of the outgoing system Y. A fourth Lorentz scalar is needed to define 
the kinematics of the hadronic part of the final state. 

The fraction of the proton momentum carried by the pomeronP] is 



where M x is the mass of the produced system. Note that the influence of t and m p 
on xp is negligible for large Q 2 and W 2 . 

The rapidity of a particle in the final state is given by 



x P = 



(p ~ p') ■ q _ M 2 x + Q 2 -t 
p ■ q W 2 + Q 2 — m 2 



(1.5) 




(1.6) 



2 We will say more about the pomeron further down. For the moment we call the object 
exchanged between the proton and the incoming probe, the pomeron. 
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Now consider the total hadronic cross section for processes like pp, pp, 7p, etc. 
as a function of the center of mass energy y/s ||. Their total cross section can be 
described as 



which is the form predicted by Regge theory. The first term describes the decrease 
of the cross section with s at low energies and the second the slow increase at high 
energies. CDF 0, Cudell et al. [|J and originally Donnachie and Landshoff || 
performed a fit using different hadronic data in each case to determine the values of 
rj and e and thus obtain a universal description of the hadronic cross section, 77 was 
found to be around 0.45, with the results for e varying on the data used for the fit 
with 0.08 from Donnachie and Landshoff to 0.112 from CDF. 

This behavior of the total cross section can be interpreted in terms of Regge 
theory M . The hadronic reaction A + B — ► C + D can be described by the t-channel 
exchange of a family of off-shell particles such that the relevant quantum numbers 
are conserved. There exists a theoretical relation between the spin J and the mass 
squared (t) for these particles which has the the approximate form 



The particles lie on so-called "Regge trajectories", with intercept aij(0) and slope a[ 
where i stands for the different particles. The dependence of the elastic cross section 
with t should behave as 



.h-h 
tot 



Ys-^ + Xs 6 . 



(1.7) 



Ji = oa{t) = «j(0) + a'jt. 



(1.8) 




(1.9) 
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The term in Eq. (|1.7|) , dominant at low energies, was fitted in Ref. || and corre- 
sponds to the intercept an ~ 0.5 for reggeon exchange, i.e. , the degenerate p, uj, f 
and a trajectories. The dominant term at high energies in Eq. (|f ,7|) corresponds to 



the so-called pomeron trajectory (see Ref. JTT]]) with an intercept ~ 1.08 according 
to the fit in ||. The slope a' P was fitted to be ~ 0.25 GeV -2 , implying that the 
exponential t distribution becomes steeper with increasing energy. 

The da/dt distribution in p - p elastic scattering has a characteristic behavior 
with an exponential fall-off, a dip and a second exponential, which is reminiscent of 
the diffractive pattern of light by a circular aperture @, hence the name Diffraction 
was used to indicate a pomeron exchange. Processes mediated by pomeron exchange 
can be classified as elastic, single diffractive and double diffractive as shown in Fig. 



L2[ From Fig. [L^, due to the lack of additional final state particles created through 
color interactions with the exchanged pomeron, it is also clear that the pomeron 
has to carry vacuum quantum numbers, i.e. , it has to be a color singlet and one 
expects to see a rapidity gapf] between the leading particle B and the system X. 

The main open questions on the pomeron are: Does one have a universal "soft" 
pomeron or a continuum of different objects with different energy dependences etc. 
up to high energies ? Can it be treated as a particle with a partonic substructure f| 
and if so what is the pomeron structure function? Does factorization hold or not? 

In the following discussion on hard diffraction and the following chapters some 
of these questions will be addressed. 



3 A region in the detector with no particle tracks. 




a) b) c) 

Figure 1.2: a) elastic scattering, b) single diffractive scattering, c) double diffractive 
scattering. 

1.2. Hard Diffraction in QCD 

Diffractive processes with one or several large scales, for example large Q 2 in deep- 
inelastic scattering (DIS), are called hard. Especially interesting are diffractive D1S 
processes at large Q 2 such as electro-production of longitudinally, polarized vector 
mesons and deeply virtual compton scattering (DVCS), which we will treat in more 
depth later, since they offer a unique opportunity to probe the interplay between 
soft and hard QCD physics. The fact that one observes a substantial number of 
diffractive like events at high Q 2 at the HERA, CDF and DO experiments, is a signal 
for the interplay of hard and soft QCD phenomena. There are two effects which 
are expected to play an important role, color transparency for systems consisting 
of quarks and gluons contained within a small size configuration and color opacity 
for large size configurations. Small size configurations lead to reactions dominated 
by hard processes with a cross section rising with energy, while reactions of large 
size configurations are dominated by soft processes. In order to test QCD on a 
quantitative and qualitative level, the ability to separate these two regimes from 
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one another becomes essential. In this context one can address the following issues: 

• The dynamics of compact systems, color transparency and perturbative color 
opacity in hard diffractive processes like electro-production of longitudinally 
polarized vector mesons or DVCS. 

• QCD predictions for the high momentum behavior of wave functions of 
hadrons consisting of light quarks and QCD physics of heavy quarkonia. 

• Violation of the DGLAP evolution equation. Due to unitarity considerations 
the increase of parton densities at small x will slow down and it is believed 
that this effect can be observed at HERA in hard diffractive processes, in the 
measurements of the proton structure functions at moderate Q 2 , as well as in 
the measurements of the structure functions of nuclei. 

• Semi-classical approximation to high energy interactions. In the limit of strong 
color fields which are typical in the HERA kinematical regime, new develop- 
ments and tests of the semi-classical approach are possible . 

1.2.1. Theoretical Foundation 

In order to understand the issue of the interplay of soft and hard physics in high 
energy reactions let us quickly classify them according to the number of scales 



involved 10 . 



Soft QCD: Soft hadron collisions are usually considered to be processes with 



a scale of ~ 1 fm. We have already discussed this topic in Sec. 1.1 



Hard QCD evolution: This second class consists of hard processes which are 
determined by two or more different scales in the interaction. DIS as well as 
hard diffractive processes belong to this class. The hard scale is provided by 
the virtuality of the photon or jets whereas the soft scale is set by the size of 
the proton. 

In order to use pQCD one has to prove that the short distance part of the 
interaction factorizes from the soft one. This is achieved for the total cross 
section and special hard diffractive processes as for example diffractive vector 



meson production and DVCS p4| , f46| by proving QCD factorization theorems. 



These factorization theorems lead to BL and DGLAP type equations |24], |46 
describing the evolution from large Q 2 to the scale of soft QCD processes, 
hence the soft physics ( in particular the soft pomeron ) enters as boundary 
conditions to these equations. 

One has to mention that the characteristic attribute of two scale processes 
is the violation of the pomeron pole factorization, i.e. , the description of 
the energy dependence of diffractive processes by a universal trajectory as 



mentioned before in Sec. [LI]. Therefore one may observe different energy 
dependences for different external particles and the energy dependence may 
change with Q 2 . 

BFKL Evolution: The original form of the BFKL equation was derived under 
the assumption of a small but fixed value of a s with the scale of a s set by the 
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scale of the external particle. The BFKL approximation therefore applies to 
processes with one large scale. HERA would offer a clean test by measuring 
high-p_|_ forward jets with p\ ~ Q 2 in low x DIS as suggested by A. Mueller 



12]. In diffractive DIS, semi-inclusive diffractive vector production with large 



momentum transfer t would be a place to look for the BFKL pomeron. 

Color Transparency and Perturbative Color Opacity: Color transparency is a 
phenomenon which describes within QCD the interaction of a small size, color 
neutral parton configuration with a hadronic target. The essence of color 



transparency is expressed by the following formula |L3 , |14| , [f0|] which follows 



from the factorization theorem for hard processes in QCD: 

of = 2a s n 2 j^Tr xG(x, 9/b 2 ), (1.10) 

where L 2 is the casimir operator of color SU(3), b is the transverse separation 
between the qq system and G stands for the gluon distribution in the target. 

The name "Color Transparency" stems from the fact that high energy pro- 
cesses are dominated by gluon exchange and that the cross section for a small 
size configuration is, through the gluon distribution, indirectly related to its 
size in the impact parameter space. The decrease of the cross section with 
decreasing size is partly compensated by the known increase of the gluon 
density. Furthermore, at fixed b, the cross section increases with increasing 
energy. Also, the interaction cross section at very high energies where there 
is a large number of gluons in the target, becomes large, naturally leading to 



perturbative color opacity. Both phenomena are of particular relevance for 
photon induced high energy interactions as the photon fluctuates into a qq 
pair and most of the time these fluctuations lead to a small size configuration 
with large k± between the pair, effectively screening each other. 

A word of caution about this classification is in order. As always, the real world 
is more complicated and therefore there is not always a clear cut distinction between 
these four types of physics. We know, for example, that the spatial size of the known 
hadrons varies from the proton with a radius of ~ 0.8 fm to the T with a radius 
of roughly 0.1 fm. Therefore, J/ty or T scattering of a proton belongs to the class 



of hard QCD evolution, as mentioned above, over a wide kinematical range fL5 |. 
However, with increasing energy the soft regime would dominate in most of the 
rapidity space of these reactions as a result of diffusion from the large scale as given 
by the mass of the heavy vector meson to the scale of soft QCD processes. 

After the above classification and comments we are now ready to supplement 
the bare statement from the previous section that diffractive events, whether hard 
or soft, are characterized by a large rapidity gap. These large rapidity gaps can 
be reasonably well described in terms of diffractive interactions as given by a phe- 
nomenological pomeron exchange. The basic idea behind this is the pomeron having 
a parton content which can be probed in hard scattering processes |§. 

For small qq configurations, which usually occur at high energies, diffractive 
scattering is driven by a two gluon coupling, which is the simplest model of a 
pomeron, since two gluons form a color neutral object. It is important to realize 

10 



that gluon radiation, which would fill the rapidity gap and coming from the pair 
of exchanged gluons, is suppressed. The issue here is that of coherence in the 
radiation from a color neutral system when the color charges are almost at the same 
impact parameter. In this case, gluon radiation with small transverse momenta, 
required to generate long range color interactions capable of filling the rapidity gap, 
is known to be suppressed |I| |TT| , since such a gluon cannot resolve a colorless 
object. Radiation of gluons with large transverse momenta is suppressed as well, 
though in this instance due to the smallness of the coupling constant. This reasoning 
is also directly applicable to hard diffractive scattering initiated by a small size qq 
pair where the exchange of a colorless pair of a hard and a relatively soft gluon can 
be calculated in QCD. These processes are still of leading twist because the QCD 



factorization theorem is modified for processes with diffractive final states fl8 , |65| . 

The above can be summarized in the following statement about the energy de- 
pendence of hard diffractive scattering: In the soft QCD regime which corresponds 
to large transverse separations of the qq pair, the parton model gives boundary 
conditions through e.g. the aligned jet model to the factorization theorem, hence 
it is substantially modified due to Q 2 evolution. In contrast to the parton model, 
one finds in QCD that the contribution of qq pairs with small b is only suppressed 
by oi s /k 2 ±_ but the interaction cross section increases rapidly at small x since it is 
proportional to the gluon density. Therefore, <t 7 *tv may increase faster with W 2 
compared to cross sections for hadron collisions since the probability of small size 
configurations in the wave function of a hadron is significantly smaller than in the 
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+ crossed diagram 




Figure 1.3: The lowest order handbag contribution to DVCS with Q 2 = —q 2 and 
q' 2 = as an example of a process where nondiagonal parton distributions are 
involved. 



photon wave function. 



1.2.2. What kind of diffractive processes are calculable in QCD ? 



Recently, it has been understood and proved ^2], |46|, ^5) that perturbative 
QCD can be applied to photon induced, inclusive, hard diffractive scattering as 



observed at HERA [^5| and photon induced, exclusive hard diffraction as for example 
electro-production of longitudinally polarized vector mesons , DVCS [[f^, [46], , 
diffractive di-muon production j|6| and high-p^, di-jet production. In chapter 4, we 
will go into the details of factorization of the DVCS and di-muon case and in chapters 
5 and 6 we will study DVCS in more detail. 

A feature particular to exclusive processes is that as a result of energy-momentum 
conservation, the fraction of proton momentum carried by the exchanged partons 
are not equal |1| |(J |26| leading to nondiagonal parton distributions (see Fig. 
OP - In order to see this more clearly think of a parton being emitted from the 
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proton with a fraction of the initial proton momentum x\ in the + direction. In 
photon induced reactions the virtual photon adds the fraction — jcjy to x\ and in the 
case of DVCS for example, the produced final state which will be measured, i.e. , 
the real photon, does not carry any + momentum away; thus the parton returning 
to the proton has momentum x\ — x which is not equal to the original momentum 
fraction x\. 

The nondiagonal distributions, their leading order evolution and their applica- 
tions particularly in the case of DVCS will be extensively studied in the remainder 
of this thesis. 

The investigation of these processes will provide novel information on the hadron 
structure and on the space-time development of QCD processes at high ener- 
gies. They offer a unique possibility to measure generalized parton distributions 
in hadrons in a new way. They will also allow us to check QCD predictions for the 
asymptotic behavior of the light-cone wave functions of hadrons. 

Chapter 2 is based on Ref. | 4"I|| , chapter 3 on Ref. | 4"0|| , chapter 4 on Ref. |4"6 



chapter 5 on Ref. chapter 6 on Ref. |TIJ and the appendix on Ref. ^2 . 
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CHAPTER 2 



Nondiagonal Parton Distributions in the Leading 
Logarithmic Approximation 



2.1. Introduction 



Due to the experimental possibility of probing nondiagonal distributions in hard 
diffractive electro-production processes, theoretical interest in this area in recent 



years |2T|, |23|, £4], |25], |42], Q has produced interesting results. A pioneering 
analysis of the nondiagonal distributions for the diffractive photo production of Z°- 
bosons in DIS where the applicability of PQCD is guaranteed was given by Bartels 



and Loewe in 1982 [27[] but went essentially unnoticed. 

In this chapter, which is heavily based on Ref . |^T| , we would like to complement 
these results by concrete predictions, albeit in the LLA, which can be tested by an 



experiment. In Sec. |27| we shall demonstrate that in the limit of small x the am- 
plitudes of hard diffractive processes can be calculated in terms of discontinuities of 
nondiagonal parton distributions. The real part of the amplitude will be calculated 
by applying a dispersion representation of the amplitude over x. We will show that 
the term in the amplitude which cannot be calculated in terms of the discontinuities 



of nondiagonal parton distributions |23|, ^4] is suppressed by one power of x in this 
limit. This result will make it possible to calculate the evolution kernels in the LLA 
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following the traditional methods |28| and to compare them to results obtained in 



the QCD-string operator approach f29 |. 



In Sec. 2.3 we calculate the nondiagonal kernels and find them equivalent to those 



in 



2^, 2IJ. They are different from the evolution equations for nondiagonal parton 



densities which were presented without derivation in |!(J. In Sec. ^| we shall make 



predictions about the nondiagonal parton distributions by solving, numerically, the 
nondiagonal GLAP evolution equations with the help of a modified version of the 
CTEQ-package. In Sec. [2.5| we shall discuss the limitations of the approximation 
and the need for NLO-results. Future directions will be discussed in the conclusions. 

2.2. Nondiagonal parton distributions and hard diffractive processes. 

It has been recently understood that the major difference in QCD between leading 
twist effects in DIS and higher-twist effects in hard diffractive processes is to be 
attributed to the fact that the latter, initiated by highly virtual, longitudinally 
polarized photons, can be calculated in terms of nondiagonal, rather than diagonal, 



parton distributions [26 



In order to calculate unambiguously hard two-body processes, it is necessary to 
calculate nondiagonal parton distributions in a nucleon. This implies knowledge of 
the non-perturbative nondiagonal parton distributions in the nucleon which have not 
been measured so far. Therefore, the aim of this section is to express the nondiagonal 
parton distributions in the nucleon through quantities being maximally close to the 
diagonal parton distributions. Our second aim is to elucidate the kinematics of the 
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nondiagonal parton distributions in the nucleon needed to describe hard diffractive 
processes. We shall also discuss the expected limiting behavior of the nondiagonal 
parton distributions. 

For the leading twist effects QCD evolution equations have traditionally been 
discussed in terms of parton distributions entering the imaginary part of the ampli- 
tude. This is because the bulk of experimental data available is on the total cross 
section of inclusive processes. This form of the evolution equation can be general- 



ized to the case of hard diffractive processes [22] and hard two-body processes |24 
The analysis of the QCD evolution equation for the nondiagonal parton densities 
shows that the evolution equation contains two terms. The first one is described by 



a GLAP-type evolution equation [22] , p3| , |24| , |2ql , whereas the second term, found in 



Ref. [Q for vector meson production at small x, cannot be interpreted in terms of 
parton distributions. The QCD evolution of this term is governed by the Brodsky- 



Lepage evolution equation |23|, E2 



2.2.1. GLAP evolution equation for hard diffractive processes 

The aim of this section is to prove that for hard diffractive processes in general, 
the (Revolution at any x in the DGLAP-region as discussed below, is described 
by a nondiagonal GLAP-type evolution equation with asymmetric D GLAP-type 
kernels and that these processes can be calculated through the discontinuity of hard 
amplitudes. This property is important for the quantitative calculations since the 
dispersive contribution has a relatively simple physical interpretation and a deep 
relation with the conventional parton densities. First, we shall deduce a relationship 
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Figure 2.1: Diagrammatic form of a factorized process in QCD. 

between amplitudes of hard two-body processes and parton densities, and we will 
find an additional term which has no probabilistic interpretation. We will restrict 
ourselves to the Q 2 -region where the parton distributions are still rising and the 
additional term is of no importance as discussed below. 

The QCD factorization theorem for hard processes (see Fig. [2.1| ) means that the 



hard part of the process can be factorized from the collinear, non-perturbative part 
up to terms suppressed by powers of Q ||34|| . The topologically dominant Feynman 
diagrams for small x processes correspond to attachments of only two gluons to the 
hard part. Although our analysis is rather general, for certainty we shall restrict 
ourselves to the case of diffractive processes where diagrams with two-gluon exchange 
dominate. [] It is convenient to decompose the momentum of the exchanged gluon 
k in Sudakov-type variables: 

k = x 1 p + pq + k t , (2.1) 



1 Hard collisions due to the exchange of 2 quarks are numerically small in the LLA at small x. 
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where 



p 2 = q 2 = and (k t q) = (hp) = 0. 



(2.2) 



To express the amplitude in terms of nondiagonal parton distributions, the contour 
of integration over (3 should be closed over the singularities of the amplitude in gluon- 
nucleon scattering at fixed X\ and x\ — x. The singularities over [3 are located in 
the complex plane of discontinuities over the virtualities of the vertical propagators: 
— and — — , and from the s- and u- channel discontinuities: ,~ %e > and j—^ — r. 

XI XI— x' (1— XI) (1+X1— X) 

The amplitude differs from depending on the different contributions of the poles 
given a certain contour of integration. This causality condition restricts the region 
of integration to: 



Our main interest is in the amplitude within the physical region where — t > but 
small as compared to other relevant scales of the process under consideration. In 
this region the amplitude can be represented as the sum of terms having s- or u- 
channel singularities only. For the s-channel contribution to the amplitude of hard 
diffractive processes, relevant for the region 1 > x± > x, the integral over f3 can only 
be closed over the discontinuity in the amplitude of the gluon-nucleon scattering in 
the variable s. 

Therefore, this contribution to the amplitude is expressed through the imaginary 
part of the amplitude for gluon-nucleon scattering. The QCD evolution of this term 
is described by a GLAP-type evolution equation where the kernel accounts for the 
off-diagonal kinematics. One also has to add a similar term corresponding to u- 



— 1 + X < X\ < 1. 



(2.3) 
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channel singularities. 

The contribution of the region x > x\ > has no direct relationship to the 
conventional parton densities. This is because the integral over (3 cannot be closed 
for s- or u-channel discontinuities but it may be closed for the discontinuities over 



the gluon "mass". In Ref. |23j the analogy of this term with the wave function of 
a vector meson has been suggested. The presence of this piece which cannot be 
evaluated in terms of parton densities introduces theoretical uncertainties into the 
treatment of hard two-body processes at large and moderate x. 

However, for the imaginary part of the amplitude, more severe restrictions on 
the region of integration apply: 

1 > Xl > x. (2.4) 

We only have to consider the discontinuity of the hard amplitude in the s- channel. 
This restriction on x\ for the s-channel contribution follows from the positivity 
requirement for the energy of the intermediate state in the s-channel cut. This 
result helps to prove that the piece which cannot be evaluated in terms of parton 
densities is inessential for hard diffractive processes. Let us now apply a dispersion 
representation over the variable s which reconstructs the real part of the amplitude 



according to the following formula for small x [33 



ReA = ImA*- dHxI ™ A) . (2.5) 
2 a m - 



X 



The only term which cannot be reconstructed in terms of a dispersion relation, 

i.e. , in terms of discontinuities of parton densities, is the subtraction constant^] in 
2 This constant is independent of s. 
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the real part. The contribution of the subtraction term to the amplitude with a 
positive signature, i.e. , symmetric under the transposition of s — > u, is suppressed 
by an additional power of s or, equivalently, by an additional power of x. For the 
processes with negative charge parity in the crossed channel, i.e. electro-production 
of a neutral pion, the amplitude is antisymmetric under the transposition?] s — > u. 
This amplitude has no subtraction terms at all, since, in QCD, it increases with 
energy slower thanf] s. Therefore, in this dispersion representation gives 

the full description. To summarize let us point out once more that the small x 
behavior of hard diffractive processes is described through the discontinuities of 
hard amplitudes. 

2.2.2. Small Xi behavior of the nondiagonal gluon distribution 

We want to stress that the slope of the x dependence of the amplitudes for diffrac- 
tive processes, however not their residue, should be independent of the asymmetry 
between fractions x\ and x<i of the nucleon momentum carried by the initial and 
final gluons. This is due to the fact that the Xi of the partons in the ladder are 
essential, but not the external x, and increase with the length of the parton ladder. 
Hence, the asymmetry between the gluons may be important in one or two rungs 
of the ladder but not in the whole ladder. Therefore, at sufficiently small x, it is 
legitimate to neglect x in most of the rungs of the ladder as compared to the x%. 

This means that the asymmetry between the gluons influences the residue but not 
3 This corresponds to a negative signature. 

4 An odderon-type contribution in PQCD is suppressed by an additional power of Q 2 . 
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the slope of the x dependence. 

Let us now discuss the small Xi behavior of g(xi,X2) - the nondiagonal gluon 



density in a nucleon. The factorization theorem - Eq. (3) of [24]- is the basis 
for the formal definition of the nondiagonal gluon density as the matrix element of 
gauge-invariant bilocal operators (cf. Eq. (6) of p4"f): 

r°° dv~~ p~ iX2P+y ~ 

g g/p (x 1 ,x 2 ,t,fi) = -/ -f - 2 (p'\TG u + (0,y-,0 T )VG» + (0)\p)(2.6) 

Here V is a path ordered exponential of a gluon field along the light-like line 
joining the two gluon operators, t is the square of the invariant momentum trans- 
ferred to the target, and /i describes the scale dependence. The sum over transverse 
gluon polarizations is implied. Actually Eq. (|2.6|) coincides with the definition given 



in 



2j| 23] for the same quantity. 
For x\ = X2, g g / p (xi,x 2 ) is related to the diagonal gluon distribution as 
xGdiag(x) = fi , nondia g (^, x). Within the leading a s \nx approximation where the 
difference between lnxj and lnx can be neglected, the nondiagonal distribution 
coincides with the diagonal one [f2l|| : 

^leading a s lnx (^1 j ^<2j t ^min' ) xCj(x'. fJ. ). (2.T) 

We want to stress here that at fixed Bjorken variable x, the cross sections of hard 
diffractive processes are expressed through g(x\, x 2 , t, /i 2 ) where x\ — x 2 = x. This 
can be proved by calculating the high energy limit of hard diffractive processes and 
then applying Ward identities similar to Ref. pi| . This means that the region of 
integration near x 2 = (x 2 <ti X\) gives only a small contribution to the amplitudes 
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Figure 2.2: Possible configuration in a parton distribution which could yield an 
enhancement for X2 — ► for a soft momentum k. 

of hard diffract ive processes. 

Alternatively, one can examine the leading regions of the integrals in the calcu- 
lation of the distribution of a parton in a parton (see Fig. [2.2| ) which is imperative in 
finding the correct hard scattering coefficients for the desired process. This calcula- 
tion is necessary since one has not only ultraviolet divergences in the partonic cross 
sections from which one wants to extract the Wilson coefficients, but also infrared 
divergences stemming from initial-state collinearities of the participating partons 



(see Ref. [|4j] for further details) which are cancelled by the perturbatively calcu- 
lated expansion of the parton distribution. The claim is that the region of X2 = 
does not give a leading contribution. This can be seen by using a simple argument 
that proper Feynman diagrams have no singularity at x<i = 0, and the region of 
integration over the exchanged gluon momenta X2 = forms an insignificant part 
of the permitted phase volume. 
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In the first step one has to show that a gluon with X2 = corresponds to a 



soft gluon and then one can use the argument by Collins and Sterman first 
introduced for proving factorization in inclusive e + e~-reactions: 

• For clarification, the quark-loop to which the gluons attach connects both the 
hard part and the part whose momenta are parallel to the vector meson and 
of the order = 

rrajv mjv 

• The minus component Z_ Q of the quark-loop is of order Q. The minus com- 
ponent of the gluon momentum k is: 



< /_. (2.8) 



~ Q 

Thus we can neglect k^ with respect to /_ in the calculation of the leading 
term of the amplitude corresponding to the leading power in the energy of the 
process. 

The transverse momentum in the quark loop (l t ) is cut off by the vector meson 
wave function and thus I 2 <C Q 2 in stark contrast to DIS, whereas the k 2 of 
the gluon is only restricted by the virtualities in the photon wave function 
which can be as high as Q 2 |. However one has to satisfy the Ward-identity 



5 We define a vector in light cone coordinates by: 

' V° + V 3 V° — V 3 



v'-<r+,v-,v l )-l-j S -,—j S -,v» 



and use the Breit frame as our frame of reference. 

6 The similarity to DIS will be restored at extremely large Q 2 as a consequence of both a 
Sudakov-type form factor in the photon vertex and a slow decrease, with increasing kt, of the 
vector meson wave function. 
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kfjT^ = where T stands for the amplitude. Using Sudakov-variables this 
becomes: 

x 2 p^ + k tp T^ = 0. (2.9) 

For x 2 — > the transverse momentum of the gluon is very small and can be 
safely neglected as compared to l t . 

• One concludes from the above said that the xi = region corresponds to a 
soft gluon {k 2 ~ 0) and we can now use the argument by Collins and Sterman, 
as will be explained in the statements below. 

Keeping the above said in mind, the /c-integral involved in the determination of 
the leading regions is of the form (up to overall factors): 

Lftk ^{{l - k) 2 + ie)(k* + te) f{l ~ k ^~ 

Lftk dk+ ((/+ - fc + )Z_ - Z? + it) (2k + k„ - k 2 + ie) f{l ~ k,P) ' (2 ' 10) 
where f(l — k,p) is the amplitude of gluon-nucleon scattering and the integrals over 

A;_ and k 2 are suppressed for convenience. If one now integrates over the remaining 

k + momentum one will have the following situations: 

• k + k_ > k 2 t : There are no obstructions in the deformation of the integration 
contour since l t ^> k t . In other words this region does not give a leading 
contribution. 

• k+k- k 2 : There are no obstructions to the contour deformation since l t 3> 
k t , hence one only has one pole in k + . Therefore, this region does not give a 
leading contribution either. 
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In conclusion, one has proved that if one of the gluons attaching the soft to the 
hard part has X2 = 0, it will be soft and thus, according to the above reasoning, 
the Xi = region of integration does not give a leading contribution to the parton 
distribution. 

2.3. Kernels in the LLA 

There are several possible ways of calculating the evolution kernels to leading order 
in QCD. We first used the traditional approach of calculating the evolution kernels 
in the LLA via the method of decay cells of e.g. a quark decaying into a quark 
and using cut-diagram techniques to calculate the appropriate Feynman graphs. 

As a cross-check we calculated the first order corrections to the bi-local quark and 
gluon operators on the light-cone, which not only yielded the nondiagonal kernels 
for the DGLAP equation but also the nondiagonal Brodsky-Lepage kernels, since we 
were calculating the whole amplitude, not only its imaginary part. However, since 
we are not interested in those kernels at the moment we will not comment on this 
fact further, let it be said though that our results on the Brodsky-Lepage kernels 
agree with those of Ref . p5[ Q . 



We performed the calculation of the cross-check in a planar gauge i.e q' ■ A — 
with q' 2 7^ 00 and used once more Sudakov-variables: 

k = pp' + aq' + k t , (2.11) 



7 Changes appropriate to the nondiagonal case were made. 

8 The advantage of such a physical gauge being that no gluons couple to the operators to first 
order, simplifying the calculations considerably. 
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where 

p' 2 = q ' 2 = and (k t q) = (k t p) = 0. (2.12) 

Since one is neglecting the proton mass one can set p' = P where P is the proton 
momentum. 

The insertion of the appropriate bi-local operators for quarks and gluons on the 
light cone into the Feynman graphs for first order corrections to those operators, 
short circuits the +-momentum in the graph, which means that the loop variable k 



has not a +-momentum of (3P but rather X\P (see |34j] for more details on calculating 
one loop corrections to parton distributions.). This fact rids us from the duty of 
taking the integral over (3. In the calculation of the kernels, it remains to take the 
integral over a which can be done by taking the residues and then isolating the 
leading term multiplying dk 2 /k 2 and the tree level amplitude. This will then yield 
the kernels in the leading logarithmic approximation.^ 

In the integral over a, one finds three different residues. Two residues stemming 
from the vertical quark or gluon propagators which yield a = ^ and a = giving 
a contribution in the Brodsky-Lepage region, i.e. , the Brodsky-Lepage kernels, and 

k' 2 

one stemming from the horizontal propagator yielding a = ^ yi J Xl ^ s which contributes 
to the DGLAP region, i.e the DGLAP kernels. This is analogous to the statements 
made in Sec. p72[ 

After having taken proper care of the definitions of our quark and gluon dis- 
tributions in the amplitudes, we find the following expressions for the nondiagonal 



9 Note that the quark to quark and gluon to gluon kernels also need the self-energy diagrams 
to regulate the two possible collinear singularities, of course, after proper renormalization. 
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evolution kernels, where A is given by x\ — x 2 and corresponds to xb^ of e.g., vector- 
meson production. For the quark — ► quark transition we find: 



v\ yi 7T 



'I- A)(l_£i) 



5(1-— )| 



(2.13) 



2/1 ' v J/1 • 

The other kernels are computed the same way from the appropriate diagrams: 



P. 
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yi vi' 



7i 
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L v{ ' j/i 



^ _|_ £1(1 _ £1^2 _ x i A l 



yi 



v{ 



(1- A)2 

v J/1 7 



[1 + (1 _ £1)2 _ A] 

q [ v j/i y j/i J 



1 



J/i 



(2.14) 
(2.15) 
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'1 _ £1^2 _l (I _ x l \ (si -A) 



J/f 7 J/l 



'1 - A)2 

J/l 7 



1 - fi 



1 _ £1 

J/i - 



zi— A 
J/i 



+ *(!--) Po 



(2.16) 



with Po = 1 1 — ^y- . A word concerning our regularization prescription is in order . In 
convoluting the above kernels, after appropriate scaling of X\ and A with yi, with 
a nondiagonal parton density, one has to replace Z\ and z 2 in the regularization 
integrals with z\ — > (yi — X\)jy\ and z 2 — > (yi — x\)/{yi — A). This leads to 
the following regularization prescription as employed in the modified version in the 
CTEQ package in the next section and in agreement with Ref. | f£^| : 

"d Vl m f 1 ^iyi/(yi)-x 1 /(x 1 ) + 

si yi yi - %\ 

1 dyt 2/1/(2/1) - x x f(x x ) f l dy x y x f{y\) - A/(ar 1 ) 



si yi 1 - X1/2/1, 

1 (m - A)/( yi ) 
«2/i- 



(2/1 -^i)(2/i - A )- 



*i 2/i 
+/(a; 1 )ln 



2/i - Zi 
1 — X\ 

1 - A 



si 2/1 



2/1- A 



(2.18) 



For A = one obtains, necessarily, the diagonal kernels, however for the distri- 
butions q = XiQ(xi, Q 2 ) and g = xiG(xi, Q 2 ), since we chose the definitions of our 
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nondiagonal distributions to go into q = xQ(x, Q 2 ) and g = xG(x, Q 2 ) rather then 



Q and g as in [|2] . We have cross-checked these results with those of Ref . [ 29 1 via 
the conversion formulas given by Radyushkin in a recent paper P5fl . The formulas 
given in a recent paper by Ji |25] do not seem to agree with ours but this is only due 
to a different choice of independent variables used by Ji. After appropriate transfor- 
mations, the formulas of E3] agree with our results ||42|| . It should be noted however 



that the kernels from Ref. |23[ are given for Q and g and not for q and g as 

we do. One just has to multiply the kernels given by those authors for the quark 
evolution equations with Xi/yi after appropriate changes for independent variables 
of course. Conversion formulas between the different notations can be found in Ref. 



The evolution equation for the quantities g(xi, A) and q(xi, A) take the following 
form in our notation: 

dg{ Xl ,A,Q 2 ) r 1 dy x r 2 2 

dq( Xl ,A,Q 2 ) f 1 dy 1 



I — \P qq q(Vu A, Ql) + P qg g( yi , A, Ql)] . (2.19) 



dlnQ 2 J xi y% 

We are interested in the calculation of the asymptotic distribution in terms of 
the symmetric distribution in the limit of small x and large Q 2 . The reason why 
this is possible is due to the fact that in this limit the main contribution originates 
from the nondiagonal distributions at X\i X2 — X\ — A with X \ 3> x\. In the case 
Xi, X2 ^> A deviations from the diagonal distribution are small and can be neglected. 

In the following section we will present the results of our numerical study and 
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show that for the case of X\ ^> x 2 — in the kinematic region of practical interest 
the diagonal and nondiagonal distribution will coincide for large Q 2 up to about a 
factor of 2. 

2.4. Predictions for nondiagonal parton distributions 

Utilizing a modified version of the CTEQ-package, we calculate the evolution of 
the nondiagonal parton distributions, starting from a low Qo = l.QGeV and with 
rather flat initial distributions for the diagonal and nondiagonal case by using the 



most recent global CTEQ-fit CTEQ4M ||. We chose g(x u x 2 ) = XiG{xi) in 



the normalization point, in accordance with our earlier argument that the possible 
difference in the distributions at small x and large Q 2 is only given by the Q 2 - 
evolution of g(xi, x 2 , t, /i 2 ). 

We have only considered light quarks, since we are interested in a proton as the 
initial state hadron and the s-quarks are only considered to give a small correction. 
The following figure (see Fig. |2.3|) shows the ratio of the nondiagonal distribution 
g(xi,x 2 ) to the diagonal distribution xG(x) from Q = 7 GeV to Q = 110 GeV and 
x 2 from ^ to x 1 with x x = 1.1 1(T 4 , 1.1 1(T 3 , 1.1 1(T 2 . 

The nondiagonal and diagonal distributions agree for x 2 — ► X\, i.e. for vanishing 
asymmetry, as expected, and within a deviation of a factor between 0.2 and 2, they 
agree for x 2 <ii. The expectation that there is no \\ix 2 contribution in the parton 
distribution, which would give a singularity for x 2 — » 0, is also supported by our 
numerical calculations. 
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Note that at large Q 2 and fixed A C 1, g(x\,X2) is determined by the initial 
parton distributions at x\, Xi ^> A where the validity of the diagonal approximation 



for gix^x-i) does not depend on our argument in Sec. 2.2. The numerical calcula- 



tion finds that the ratio of nondiagonal to diagonal distribution is larger than 1 as 



anticipated by Radyushkin |38[ based on general arguments about the nature of the 
double distribution which he discusses in pT2"|j . 

To see whether our numbers, i.e. , our numerical methods, could be trusted, 
we used a MATHEMATICA program to calculate the first iteration and the first 
derivative of the evolution to see how good or bad our numbers were. As it turns 
out our integration routines produce a very good agreement with the numbers from 
MATHEMATICA with a relative difference of 5%. This leads us to believe that our 
numbers can be trusted to high accuracy for x 2 of 0[x,\) and within 5% at X2 down 
by two orders of magnitude as compared to X\. 

A few words about the nature of the modifications to the CTEQ-package are in 
order at this point. The basic idea we employed, was the following: In the CTEQ 
package the parton distributions are given on a dynamical x- and Q-grid of variable 
size where the convolution of the kernels with the initial distribution is performed on 
the x-grid. Due to the possibility of singular behavior of the integrands, we perform 
the convolution integrals by first splitting up the region of integration according to 
the number of grid points, analytically integrating between two grid points Xi and 
Xi + \ and then adding up the contributions from the small intervals. We can do 
the integration analytically between two neighboring grid points by approximating 
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the distribution function through a second order polynomial ay 2 + by + c, using the 
fact that we know the function on the grid points Xi-\,Xi and Xi + \ and can thus 
compute the coefficients a,b,c of the polynomial. This approximation is warranted 
if the function is well behaved and the neighboring grid points are close together. 
We treat the last integration between the points Xi and x<i (which are not to be 
confused with the x\ and X2 of the parton ladder) by taking the average of x\ and 
X2 and the values of the function at x\ and X2 and using those averages together 
with X\, X2 and the value of the function at X\ and X2 to compute the coefficients 
of the polynomial^- The coefficients are computed in the new subroutine NEWAR- 
RAY and the integration of the different terms in the kernels is performed in the new 
subroutine NINTEGR. Appropriate changes in the subroutines NSRHSM, NSRHSP 
and SNRHS were made to accommodate the fact that the kernels and also the inte- 
gration routines changed from the original CTEQ package. A detailed description 
of the code will be in the next chapter; also see Ref. [HDJ . 

2.5. Limitations of the LLA in the nondiagonal case 

The LLA approach of the previous sections accounts for the contribution of a certain 
rather limited range of integration in the parton distributions. Regions outside 
these limits might contribute to the leading power. Looking at some other physical 
quantities such as F 2 , where one finds substantial modifications due to the NLO- 

terms, we are forced to assume that this may be also true in our case. This results in 

10 See the next chapter for an updated prescription. The results of the code are, however 
unchanged. 
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the urgent need to carry out a NLO calculation and numerical study of the evolution 
equation, which will be the next step of our program. 

2.6. Conclusions and Outlook 

In summary, we have calculated the evolution kernels for nondiagonal parton dis- 
tributions in the LLA using traditional methods and found agreement with the 



results of |23|, |25], ^9f deduced by other methods. It was important to show that 
the traditional approach can still be applied and thus traditional methods can be 
used to calculate systematically hard diffractive processes within the NLO approxi- 
mations. We have also proved the similarity between the diagonal and nondiagonal 
parton distributions. The latter ones determine the cross sections of hard diffractive 
processes in the small x region. We have made predictions about the nondiagonal 
parton distributions within the LLA with the help of a modified version of the 
CTEQ-package. Numerical calculations found the diagonal and nondiagonal gluon 
distributions, which dominate hard diffractive processes, to be very similar at small 
x as expected from the previous discussion. 
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Figure 2.3: The fraction g{x\,X2)/x-\G(x\) as a function of £2 for fixed 
x\ and Q = 9.3 (solid line), 17.1 (dotted line), 40.3 (dashed line) and 
110 (dashed-dotted line) GeV 2 . 
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CHAPTER 3 



Methods in the LO Evolution of Nondiagonal Parton 
Distributions: The DGLAP Case 



3.1. Introduction 



In this chapter, which is based on Ref. [fEJ, we give an exposition of the algorithms 



used to numerically solve the generalized GLAP-evolution equations of the last 
chapter. The main part of the evolution program was taken over from the CTEQ 
package for the diagonal parton distributions from inclusive reactions. At this point 
in time the evolution kernels for generalized parton distributions are known only to 
leading order in a s , as pointed out previously and thus our analysis will be a leading 
order one. 

This chapter is organized in the following way. In Sec. |3.2| we will quickly re- 
view the formal expressions for the parton distributions and the evolution equations 
together with the explicit expressions for the kernels and a first comment on the 
resulting numerical problems. In Sec. |3.3| we will explain the difference between 
our algorithms and the ones used in the original CTEQ package and then give a 
detailed account of how we implemented our algorithms. In Sec. |3.4| we demonstrate 
the stability of our code and show that we reproduce the case of the usual or di- 



agonal parton distributions within 0.5% for a vanishing asymmetry factor. Sec. |3T5 
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+ crossed diagram 




Figure 3.1: The lowest order handbag contribution to DVCS with Q 2 = —q 2 and 
q' 2 = 0. 

contains concluding remarks. 

3.2. Review of Nondiagonal Parton Distributions, Evolution Equations 
and Kernels 

3.2.1. Nondiagonal Parton Distributions 

As explained in the previous chapter, generalized or nondiagonal parton distribu- 
tions occur for example in exclusive, hard diffractive J/ip or p meson production or 
alternatively in deeply virtual Compton scattering (DVCS), where a real photon is 
produced^. As mentioned in Sec. 3T, since one imposes the condition of exclusive- 



ness on top of the diffraction condition, one has a kinematic situation in which there 
is a non-zero momentum transfer onto the target proton as evidenced for example 



by the lowest order "handbag" diagram of DVCS in Fig. |3J . 

The nondiagonal quark and gluon distributions have the following formal def- 
inition as matrix elements of bilocal, path-ordered and renormalized quark and 



'We will say more about DVCS in the following chapters. See also Ref. g|, g| gg, |1J |§ |§ 

[ n n M 
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gluon operators sandwiched between different momentum states of the proton as 
in the factorization theorems for exclusive vector meson production pi| and DVCS 

WA M 



U/p= r ^e-^ + ^(p|T^(0, 2/ -,O ± ) 7 + ^(0)b'), 

J-oo 47T 

h/v = - [°° T ^e-^ + ^( J 9|TG+(0, Z/ -,0 ± )PG^(0)b'). (3.1) 

J —oo ATT X\%2P 

with X2 = x\ — A where the asymmetry or nondiagonality parameter A is usu- 
ally XBj in, for example, DVCS or exclusive vector meson production however not 
in diffract ive di-muon production where A = xgj + £1 and £1 is the longitudinal 
momentum fraction of the produced time-like photon decaying into a yU + yU~-pair. 

3.2.2. The GLAP-Evolution Equations for Nondiagonal Parton Distribu- 
tions 

Reprising on the discussion of chapter 2, the GLAP-evolution equations follow from 
the usual renormalization group transformation of the parton distributions and lead 
to the following evolution equations for the singlet (S) and non-singlet (NS) case 



41, 25, 42 



dq NS (x u A,Q 2 ) 


f 1 dm p 


dlnQ 2 


Jxi yi 


dg s (xi, A,g 2 ) 


r 1 dm - 


d\nQ 2 


Jxi yi L 


dq s (xi, A,g 2 ) 


f 1 dyi - 


d\nQ 2 





ggtto\tiL, — , -su/ 1 gq^ 



P q b q q s ( yi , A > Qo) + PaMVh a, Qi) ■ (3.2 



Note that qs,NS — x ifq/p, 9s — x i x 2f g / p and the kernels to leading order] 2 ] were 

given in the previous chapter together with a discussion on the the generalized + 
2 For more details on the derivation of the kernels to leading order see for example ^||, f42f . 
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regularization prescription. 

3.3. Differences between the CTEQ and our Algorithms 

Let us point out in the beginning that our code is to 99% the original CTEQ-code 
(for more on this code see Ref. |4l|). We only modified the subroutines NSRHSM, 
NSRHSP and SNRHS within the subroutine EVOLVE and added the subroutines 
NEWARRAY and NINTEGR. These routines are only dealing with the convolu- 
tion integrals but not with, for example, the (^-integration or any other part of 
the CTEQ-code which remains unchanged. This is due to the fact that the main 
difference between the diagonal and nondiagonal evolution stems from the different 
kernels which only influence the convolution integration and nothing else. 

In order to make the simple changes in the existing routines more obvious we 
will first deal with the new subroutines. 

3.3.1. NEWARRAY and NINTEGR 

Due to the increased complexity of the convolution integrals as compared to the 
diagonal pointed out in Sec. [3.2.2| , we were forced to slightly change the very 



elegant and fast integration routines employed in the original CTEQ-code. The basic 
idea, very close to the one in the CTEQ-code, is the following: Within the CTEQ 
package, the parton distributions are given on a dynamical x- and Q-grid of variable 
size where the convolution of the kernels with the initial distribution is performed on 
the x-grid. Due to the possibility of singular behavior of the integrands, we perform 
the convolution integrals by first splitting up the region of integration according 
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to the number of grid points in x, analytically integrating between two grid points 
Xi and Xi + \ where i runs from 1 to the specified number of points in x and then 
adding up the contributions from the small intervals as exemplified in the following 
equation: 

f 1 d ^Kf{x 1 /y 1 ^/y 1 ,y 1 ) = Sf =0 T +1 —Kf(x 1 /y 1 ,A/y 1 ,y 1 ), (3.3) 
J xi yi Jx % yi 

where Kf(xi/yi, A/yi, yi) is the product of the initial distribution for each evolution 
step and an evolution kernel with x = %n = 1- We can do the integration 
analytically between two neighboring grid points by approximating the distribution 
function f(yi) through a second order polynomial ayf + by% + c, using the fact that 
we know the function on the grid points x« and x i+ i and can thus compute the 
coefficients a,b,c of the polynomial in the following way, given the function is well 
behaved and the neighboring grid points are close together Q 

f(x 1+ x) = ax 2 i+1 + bx i+1 + c 

f(xi) = ax\ + bxi + c 
fixi-t) = axl_ x + + c (3.4) 

which yields a 3 x 3 matrix relating the coefficients of the polynomial to the values 
of the distribution functions at Xi-\,Xi and a^+i. Inverting this matrix in the usual 
way one obtains a matrix relating the x values of the distribution function to the 

coefficients making it possible to compute them just from the knowledge of the 

3 The parton distributions functions are smooth and well behaved thus one just has to use 
enough points in x. 
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different x values and the value of the distribution function at those x values. This 
calculation is implemented in NEWARRAY where the initial distribution is handed 
to the subroutine and the coefficient array is then returned. The coefficient array in 
which the values of the coefficients for the integration are stored, has 3 times the size 
of the user-specified number of points in x since we have 3 coefficients for each bin 
in x. We treat the last integration between the points x$ and X\ by approximating 
the distribution in this last bin through a second order polynomial. However, for 
this last bin, the coefficients are computed using the last three values in x and of the 
distribution at those points, since the point x_\ which would be required according 
to the above prescription for calculating the coefficients, does not exist. 

After having regrouped the terms appearing in the convolution integral in such 
a way that all the necessary cancelations of large terms occur within the analytic 
expression for the integral and not between different parts of the convolution integral, 
the integration of the different terms is performed in the new subroutine NINTEGR 
with the aid of the coefficient array from NEWARRAY. 

As mentioned above the convolution integral from x\ to 1 is split up into several 
intervals in which the integration is carried out analytically. To give an example 
of this procedure we consider the convolution integral of P qg (x\/y\, A/yi) with the 
parton distribution gs{lJi)'- 

f 1 dyx r 1 dyi x\{xi - A) r l dy l xt(yi - x x ) 2 

J*i yi J ">i yi yi{yi-A) J ^ yi 2/1(2/1 — ^) 

suppressing presently irrelevant factors in front of the integral. The two parts in Eq. 



(|3.5| ) are calculated in different parts of NINTEGR and then put together in either 
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NSRHSM, NSRHSP or SNRHS. 

In NINTEGR the integrals are split up according to Eq. (|3.3| ) and then ana- 
lytically evaluated in the different x-bins []. If the dependence of the integrand on 
A is only of a multiplicative nature it is enough to compute the integral for each 
bin once. To get the value of the convolution integral for a term with such a A | 
dependence, it is enough to store the result of the integration in the bin from xn~i 
to xn in the output array for this term at the position iV — 1 add to this result the 
value of the integral in the bin from to x/v-i and store it at the position N — 2 
and so forth. In this manner one only has to calculate N — 1 integrals, however 
if the integrand has a more complicated dependence on A like x± — A one needs 
to compute N(N — l)/2 integrals. For example in order to find the integration 
value for the x^-i bin with x\ = x^-i one needs only one integral but at x^-2 we 
have to redo our integral for the x^-i bin since x\ = x^-2 plus we need to add 
the contribution from the xn-2 bin to get the correct answer for the output array 
at position N — 2 and so forth. This need for additional evaluations of integrals 
slows the program down but in the end it turns out to be only about a factor of 
4 — 5 slower than the original CTEQ-code which is speed optimized. The integral 
with the regular + - prescription is evaluated using the routine HINTEG from the 
original CTEQ-code whereas the generalized + - prescription is evaluated according 

to the methods described above due to its nontrivial dependence on x\ and A. 

4 The general analytic expressions for the convolution integrals in an arbitrary cc-bin were ob- 
tained with the help of MATHEMATICA. 
5 The value of A is specified in NINTEGR. 

6 The value of the output at position N is always since in this case the upper and the lower 
bound of the integral coincide. 

40 



In the case of A << x\ the analytic expressions obtained for the above general 
case are expanded to first order in A and then the same methods as above for 
evaluating the integrals are applied. The last case also allows us to go to the 
diagonal case by setting A = without using the integration routines from the 
original CTEQ-code giving us a valuable tool to compare our code to the original 
one. 

3.3.2. Modifications in NSRHSM, NSRHSP and SNRHS 

The modification in the already existing routines NSRHSM, NSRHSP and SNRHS 
of the original CTEQ package are rather trivial. The most notable difference is that 
the subroutine NEWARRAY is called every time either of the three subroutines is 
called since the distribution function handed down on an array changes with every 
call of NSRHSM, NSRHSP and SNRHS. In NSRHSM and NSRHSP, NEWARRAY 
is only called once since one is only dealing with the non-singlet part containing no 
gluons, whereas in SNRHS the subroutine for the singlet case, one needs a coefficient 
array for both the quark and the gluon. Besides this change, the calls for INTEGR 
are replaced by NINTEGR according to how the convolution integral has been 
regrouped as explained in Sec. |3.3.1|. The different regrouped expressions are then 



added, after integration for different x- values, to obtain the final answer in an output 
array which is handed back to the subroutine EVOLVE. The method is the same as 
in the original CTEQ-code but the terms themselves have changed of course. 



41 



3.4. Code Analysis 

As a first step we tested the stability and speed of convergence of the code and 
found that by doubling the number of points in the x-grid, which is only relevant 
for the convolution integral, from 50 to 300 the result of our calculation changed by 
less than 0.5%, hence we can assume that our code converges rather rapidly. We 
also found the code to be stable down to an x 2 = 10~ 10 beyond which we did not 
test. Furthermore we can reproduce the result of the original CTEQ-code, i.e. the 
diagonal case in LO within 0.4% giving us confidence that our code works well since 
the analytic expressions for the diagonal case are the expansions of the general case 
of non- vanishing asymmetry up to, but not including, 0(A 2 ). 

In the following figures (Fig. |3~^ - |3.7| ) we compare, for illustrative purposes, the 



diagonal and nondiagonal case by plotting the ratio 

g(xi,x 2 ,Q 2 



R g (x 1 ,x 2 ,Q d 



x 1 G(x 1 , Q 



2^ 



2 q(xi 1 x 2 ,Q 2 ^ 



for various values of x±, Q 2 and A = xsj [Ji-e. , varying x 2 , using the CTEQ4M and 
CTEQ4LQ f\ parameterizations [p9 |. We assume the same initial conditions for the 



diagonal and nondiagonal case (see Ref. ||T for a detailed physical motivation of 
this ansatz). 

The reader might wonder why only CTEQ4M and CTEQ4LQ and not GRV 



7 We also plot the same ratio for A = to demonstrate the deviation from our code in the 
diagonal limit from the CTEQ-code. 

8 CTEQ4LQ gives the best fit at low Q 2 whereas CTEQ4M gives the best x 2 -ht for a large 
range of Q and x. 

42 



or MRS were used. The answer is not a prejudice of we against GRV or MRS but 
rather the fact that a comparison of CTEQ4M and CTEQ4LQ shows the same char- 
acteristic as comparing, for example, CTEQ4M and GRV at LO. The observation 
is the following: CTEQ4LQ is given at a different, rather low, Q, as compared to 
CTEQ4M and hence one has significant corrections from NLO terms in the evolution 
at these scales. This leads to a large difference between CTEQ4LQ and CTEQ4M 



(see Fig. |3T8| ), if one evolves the CTEQ4LQ set from its very low Q scale to the 
scale at which the CTEQ4M distribution is given, making a sensitivity study of 
nondiagonal parton distributions for different initial distributions impossible at LO. 
Of course, the inclusion of the NLO terms corrects this difference in the diagonal 
case but since there is no NLO calculation of the nondiagonal case available yet, a 
study of the sensitivity of nondiagonal evolution to different initial distributions has 
to wait. 

The figures themselves suggest the following. First, the lower the starting scale, 
the stronger the effect of the difference of the nondiagonal evolution as compared to 
the diagonal one and also that most of the difference between nondiagonal and diago- 
nal evolution stems from the first few steps in the evolution at lower scales. Secondly, 
under the assumption that the NLO evolution in the nondiagonal case will yield the 
same results for the parton distributions at some scale Q, irrespective of the starting 
scale Q , in analogy to the diagonal case. One can say that the NLO corrections to 
the nondiagonal evolution will be in the same direction and same order of magnitude 
as the diagonal NLO evolution. If, in the nondiagonal case, the NLO corrections 
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were in the opposite direction, which would lead to a marked deviation from the 
LO results, compared to the diagonal case, the overall sign of the NLO nondiagonal 
kernels would have to change for some A / since in the limit A — > we have 
to recover the diagonal case. This occurance is not likely for the following reason: 
First, the Feynaman diagrams involved in the calculation of the NLO nondiago- 
nal kernels are the same as in the diagonal case, except for the different kinematics, 
therefore, we have a very good idea about the type of terms appearing in the kernels, 
namely polynomials, logs and terms in need of regularization such as ln(z) hi ^~*) • 
Moreover, the kernels, as stated before, have to reduce to the diagonal case in the 
limit of vanishing A which fixes the sign of most terms in the kernel, thus the only 
type of terms which are allowed and could change the overall sign of the kernel are 
of the form 

-f(x 1 /y 1 ,A/y 1 ) (3.7) 
Vi 

which will be numerically small unless y% ~ A in the convolution integral of the 
evolution equations. Moreover, we know that in this limit the contribution of the 
regularized terms in the kernel give the largest contributions in the convolution 
integral and therefore sign changing contributions in the nondiagonal case would 
have to originate from regularized terms. This in turn disallows a term like Eq. 
(|3.7|) due to the fact that regularized terms are not allowed to vanish in the diagonal 
limit, since the regularized terms arise from the same Feynman diagrams in the both 
diagonal and nondiagonal case. Therefore, the overall sign of the contribution of 
the NLO nondiagonal kernels will be the same as in the diagonal case. 
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A word should be said about how the results of Ref. [PJ] compare to ours. For 
the case of the same A = 10 -3 similar starting scales and almost identical values of 
Q we find good agreement with their numbers for R g at X\ ~ A [] and are slightly 
higher at larger x\. The observed differences are due to the fact that the quark 



distributions are included in our evolution as compared to MM and their initial 



distributions are slightly different. We also find very similar ratios to MM if one 



changes the starting scale to a lower one. The slight difference of a few percent in 



the ratios between us and J44J can again be attributed to the fact that they used 
the GRV distribution as compared to our use of the CTEQ4 distributions, hence 
a slight difference in the starting scales and their lack of incorporating quarks into 
the evolution. 

3.5. Conclusions 

We modified the original CTEQ-code in such a way that we can now compute the 
evolution of nondiagonal parton distributions to LO. We gave a detailed account of 
the modifications and the methods employed in the new or modified subroutines. 
As the reader can see, the modifications and methods themselves are not something 
magical but rather a straightforward application of well known numerical methods. 
We further demonstrated the rapid convergence and stability of our code. In the 
limit of vanishing asymmetry we reproduce the diagonal case in LO as obtained 
from the original CTEQ-code within 0.4%. We also have good agreement with the 



9 This was also the case in Ref. where we initially put the energies as Q 2 where in fact they 
are given as Q, which led to some confusion in the comparisons of this first study to Ref. EH. 
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results in Ref. [44|. In the future, after the NLO kernels for the nondiagonal case 
have been calculated, we will extend the code to the NLO level to be on par with 
the diagonal case. 
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Figure 3.2: R g is plotted versus x\ for fixed A using the CTEQ4M parameterization 
with Qq =1.6 GeV and A = 202 MeV. 
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Figure 3.3: R q is plotted versus x\ for fixed A using the CTEQ4M parameterization 
with Qq =1.6 GeV and A =202 MeV. 
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Figure 3.4: R g and R q are plotted versus x\ for A = using the CTEQ4M param- 
eterization with Q =1.6 GeV and A = 202 MeV. 
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Figure 3.5: R g is plotted versus x\ for fixed A using the CTEQ4LQ parameterization 
with Q =0.7 GeV and A =17 A MeV. 
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Figure 3.6: R q is plotted versus x± for fixed A using the CTEQ4LQ parameterization 
with Q =0.7 GeV and A =17 A MeV. 
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Figure 3.7: R g and R q are plotted versus x\ for A = using the CTEQ4LQ param- 
eterization with Q =0.7 GeV and A =174 MeV. 
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Figure 3.8: The ratios for CTEQ4M to CTEQ4LQ for gluons and quarks in the 
diagonal case is plotted to demonstrate the difference between the LO evolution for 
these parameterizations. 
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CHAPTER 4 



Proof of Factorization for Deeply Virtual Compton 

Scattering in QCD 



4.1. Introduction 



In this chapter, based on Ref. J46|, we prove factorization for the deeply virtual 
Compton scattering (DVCS) amplitude in QCD up to power suppressed terms, to 
all orders in perturbation theory. This proof is important because of the recent great 



interest in DVCS f2|, g|, |§ [5T|, ||, |53], |4], |5§ |§. One important use of DVCS is 
as a probe of off- forward (or nondiagonal) distributions [25, 49, 57, 41 1. These 
differ from the usual parton distributions probed in inclusive reactions by having a 
non-zero momentum transfer between the proton in the initial and final state. 
A related process which is also used to probe off-diagonal parton densities is 



exclusive meson production in deep-inelastic scattering [21| pq| , for which a proof 



of factorization was given in ||24j| . Compared with this process, DVCS is simpler 
because the composite meson in the final state is replaced by an elementary particle, 
the photon, and thus there is no meson wave function in the factorization formula. 

Ji [^5j and Radyushkin j|9| have provided the key insights that indicate that a 
factorization theorem is valid for DVCS and Radyushkin also gave an all-orders proof 



49| . In this chapter, we provide an alternative proof, and give a new treatment of 
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some problems that were touched upon in Ref. |49| but that were not fully solved. 
The proof follows the general lines of proofs of factorization for other processes 
|H. 35], and the most noteworthy feature is that the proof is simpler than for any 



other process. Even for ordinary deep-inelastic scattering one needs to discuss the 
cancelation of soft gluon exchanges and of final-state interactions, whereas these 
complications are not present in the leading power for DVCS. 

The chapter is organized in the following way: First we will state the theorem 



to be proved in Sec. £12| and then explain the steps necessary to prove it, as adapted 



from Ref. f24j| , in Sec. |4.3| . The complications mentioned above concern the situation 



when one of the two lines connecting the parton density to the hard scattering carries 



zero longitudinal momentum, and these are given a detailed treatment in Sec. [4.3.8 . 
In the last section we will make concluding remarks. 

4.2. Factorization Theorem 

The process under consideration is DVCS, which is the elastic scattering of virtual 
photons: 

-f*(q)+P(p)^l { *\q')+P'(p-A) (4.1) 

where the diffracted proton P' may also be replaced by a low-mass excited state and 
the final-state photon can be either real or time-like. This process is the hadronic 
part of ep — *> e^p' for a real photon or of ep — > efi + fi~p f for a time-like photon. 
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It is convenient to use light-cone coordinates with respect to the collision axisQ. 
The momenta in the process then take the form: 



2 

m 



q - \ P '2xp+' '-)• 

^ ( J ( 1 + ^- 2(l-x- M )^ -^J- (4 - 2) 
Here, x is the Bjorken scaling variable, Q 2 is the virtuality of the initial photon, m 2 
is the proton mass, t = A 2 is the momentum transfer squared, and a is a parameter 
that specifies the virtuality of the outgoing photon: q' 2 = aQ 2 . Thus, a = for 
a real photon and a > for a time-like photon. Finally, ~ means "equality up to 
power suppressed terms". 

The theorem to be proved is that the amplitude for the process (|4.1|) is: 

T = Ej / dxifi/p {x\, x 2 , t, /i) Hi (xi/x, x 2 /x, //) + power-suppressed corrections, 

J — l+x 

(4.3) 

where the fi/ p is a nondiagonal parton distribution and Hi is the hard-scattering 

coefficient for scattering off a parton of type i. We let X\ be the momentum fraction 

of parton i coming from the proton, so that x 2 = X\ — (1 + oc)x is the momentum 

fraction which is returned to the proton by the other parton line joining the parton 
1 We define a vector in light cone coordinates by: 

/ V° + V 3 V° — V 3 
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distribution and the hard scattering. There is implicit polarization dependence in 
the amplitude. /i is the usual renormalization/factorization scale which should be 
of order Q to allow calculations of the hard scattering coefficients within finite-order 
perturbation theory. The \x dependence of fi/ p is given by equations of the DGLAP 
and Brodsky-Lepage kind 09l [57], Ell |24|, HE. The parton distributions in Eq. 



together with their evolution equations, are defined using the conventions of 
p4] , PH . They may easily be transformed into those given in [ 25 , Ej| by a change of 



normalization and of kinematic variables. 
4.3. Proof of Theorem 

The proof of our theorem Eq. ( |4.3| ) can be summarized as follows |: 

• Establish the non-ultra-violet regions in the space of loop momenta contribut- 
ing to the amplitude. 

• Establish and prove a power counting formula for these regions. 

• Determine the leading regions of the amplitude. 

• Define the necessary subtractions in the amplitude to avoid double counting. 

• Taylor expand the amplitude to obtain a factorized form. 

• Show that the part containing the long-distance information can be expressed 
through matrix elements of renormalized, bi-local, gauge invariant operators 
of twist-2. 



For a very detailed account of the basic steps and potential problems see Ref. 
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4.3.1. Regions 

First let us establish the regions in the space of loop momenta contributing to the 
asymptotics of the amplitude, i.e., the generalized reduced graphs. The steps leading 
to the generalized reduced graphs are identical to the steps 1-3 in Sec. IV of Ref. 



24| , i.e., scale all momenta by a factor Q/m, use the Coleman-Norton theorem to 
locate all pinch-singular surfaces in the space of loop momenta (in the zero-mass 
limit), and finally identify the relevant regions of integration as neighborhoods of 
these pinch-singular surfaces. 



In the first step, the scaling of momenta, we proceed as follows 59 . We write 



the general momentum and a general mass m in units of the large scale Q: 

P = Qjfe", m = Qm. (4.4) 

Due to working in the rest frame of the virtual photon both of the light like compo- 
nents are of order Q. Therefore, when everything is expressed in terms of the above 
scaled variables, dimensional analysis shows that the large Q limit is equivalent to 
the rh — > limit. Since the amplitude is dimensionless we have 

T(Q 2 ,p, q', A, m; y) = T{l,p, q', A, rh, Q/y), (4.5) 

by regular dimensional analysis. 

The most basic region is found where all internal lines obey k 2 > Q 2 , with 
the scaled momenta k having virtualities of 1 or bigger. In such a region one is 
entitled to setting the masses equal to zero, make the external hadrons light-like 
and set the renormalization scale /i equal to Q, thus avoiding large logarithms. As 
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it turns out, however, this region is not only not the only one but it is even not 
leading. Nevertheless, one can now see that all other relevant regions correspond to 
singularities of massless Feynman graphs. They are neighborhoods of pinch-singular 
surfaces of massless graphs, in other words, surfaces where the loop momenta are 
trapped at singularities. The conditions for a pinch singularity are the Landau 
conditions for a singularity of a graph. Only pinch singularities are of interest, since 
at a non-pinched singularity, one can deform the contour of integration such that at 
least one of the singular propagators is no longer near its pole. In case of a pinch 
singularity from propagator poles in the massless limit, we know then that in the 
real graph, with nonzero masses but large Q, the contour of integration is forced 
to pass near the propagator poles. Consequently, it is not possible to neglect the 
masses in this region. Conversely, if the contour is not trapped by the poles, the 
contour can be deformed away from the poles, and the mass may be neglected in 
evaluating the corresponding propagators. 

In the second step we use the Coleman-Norton theorem [6(J which states that 
each point on a pinch- singular surface in loop momentum space, corresponds to 
a space-time diagram obtained in the following way. First one obtains a reduced 
graph by contracting to points all of the lines whose denominators are not pinched. 
Then one assigns space-time points to each vertex of the reduced graph in such a 
way that the pinched lines correspond to classical particles. In other words each line 
is assigned a particle propagating between space-time points corresponding to the 
vertices at its ends. The momentum of the particle is exactly the momentum carried 
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by the line, with its orientation such that it has positive energy. If, for some set 
of momenta, one cannot construct such a reduced graph, one is free to deform the 
contour of integration. A reduced diagram, therefore, corresponds to a classically 
allowed space-time scattering process. In the zero mass limit, the construction of 
reduced graphs becomes very simple, since all pinched lines must carry either light- 
like or zero momentum. Furthermore, each light-like momentum must be parallel 
to one of the, now light-like, external lines. 

To be precise, in the zero mass limit of the process under consideration we have: 

• One light-like incoming proton of momentum p\ = (p + , 0, 0j_). 

• One light-like outgoing proton line of parallel momentum 
p'Z = ((l-{l + a)x)p+,Q,0 ± ). 

• One light-like outgoing photon line of momentum q ffJ- = (axp + , Q 2 /2xp + , 0j_). 

• One incoming virtual photon of momentum = (—xp + , Q 2 /2xp + , 0±). 

The results of the above construction are the two kinds of reduced graph shown 



in Fig. O . There, A and B denote collinear graphs with one large momentum 
component in the + and — direction respectively, H denotes the hard scattering 
graph, and S denotes a graph with all of its lines soft, i.e., in the center-of-mass 
frame all the components of the momenta in S are much smaller than Q. Note that, 
of the external momenta, p and p' belong to A, q' belongs to B or H , and q belongs 
to H. 
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a) b) 

Figure 4.1: a) Reduced graph for DVCS with direct coupling for the out-going 
photon to hard subgraph, b) The same without a direct coupling for the out-going 
photon. 

When the two external photons have comparably large virtualities, the only 



reduced graphs are of the first kind, Fig. [lT| a, where the out-going photon couples 



directly to the hard scattering. But when the outgoing photon has much lower 
virtuality than the incoming photon, for example, when it is real, we can also have 
the second kind of reduced graph, Fig. P~l~|b, where the out-going photon couples to 



a B subgraph. As we will see later, power counting will show that the second kind 
of reduced graph, Fig. P~T|b, is power suppressed compared to the first kind, with a 



direct photon coupling. This implies that we will avoid all the complications which 
were encountered in P41] that are associated with the meson wave function. 



The corresponding space-time diagram is Fig. [4.2| . In this figure, each solid line 
corresponds to a light-like line of the reduced graphs, with a 45° orientation to 
correspond to their light-like lines of propagation. The dashed lines correspond to 
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a) b) 

Figure 4.2: a) Space-time picture of the DVCS process with a collinear-to--B part 
as in vector meson production, b) Space-time picture of the DVCS process without 
a collinear-to--B part. 

the soft subgraph S. As far as the Coleman-Norton theorem is concerned the lines 
are degenerate. In fact they are carrying zero-momentum implying that they have 
no specific orientation. They are therefore indicated by curved lines of no particular 
orientation. The location of the endpoints of the soft lines can be anywhere along 
the world lines of the collinear lines. The hard vertex H occurs at the intersection of 
the collinear lines. The world line , in the + direction, of the collinear-to-A subgraph 
actually consists of several lines propagating together and possibly interacting with 
each other. 

In the space-time representation of a Feynman graph, there is normally an ex- 
ponential suppression when there are large space-time separations between vertices. 
One obtains a singularity when this suppression fails and the Coleman-Norton con- 
struction gives exactly the relevant configurations of the vertices. The singularity is 
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generated by the possibility of integrating over arbitrarily large scalings^ in coordi- 
nate space without obtaining exponential suppressions. 

Note that the above discussion relies on the use of a covariant gauge. The 
use of an axial gauge, albeit convenient, leads to unphysical singularities in the 
propagators. These singularities do not give the normal rules of causal relativistic 
propagation of particles and, furthermore, make the derivation of the factorization 



theorem, beyond leading log, very difficult |35|, |36 
4.3.2. Power Counting 

Each reduced graph codes a region of loop-momentum space, a neighborhood of 
the surface 7r of a pinch singularity in the zero-mass limit. The contribution to the 
amplitude from a neighborhood of n behaves like Q p ^\ modulo logarithms, in the 
large-Q limit, with the power given by 

p(7r) = 4 — n(H) — ^(quarks from S to A, B) — 3#(quarks from S to H) 

-2#(gluons from S to H). (4.6) 

where n(H) is the number of collinear quarks, transversely polarized gluons, and 
external photons attaching to the hard subgraph H. Such results were obtained 
by Libby and Sterman [[5^, |6l| . The particular form of Eq. (|4.6|) was given in |24 



together with a proof that applies without change to DVCS. 

We will detail it here, nevertheless, for completeness sake. The arguments used 
in the proof will rely on general arguments about dimensional analysis and Lorentz 



3 The scaling of the world lines in a reduced graph by a common factor does not affect the 
properties of that graph. 
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boosts. 

We first consider the case of only the hard and collinear subgraphs without a 
soft subgraph. Let the hard subgraph have N q external quark/antiquark lines and 
N g external gluons, as well as two photon lines. By definition, all components in the 
hard subgraph have virtuality of order Q 2 . Since the hard subgraph has dimensions 
du = 2 — ~N q — N g and all the couplings are dimensionless, it contributes a power 

Qd„ = Q2-lN q -N 9 ( 47) 

to the amplitude^. 

For the momenta collinear to the proton we have 

(TV? \ 
Q, ml . (4-8) 

Since x is small, there are also collinear momenta with + components much larger 
then Q. We will deal with this problem later on; let us just assume for the moment 
that x is not small. 

The collinear configurations can be obtained by boosts from a frame in which all 
components of all momenta are of order m. Since the virtualities and the sizes of 
regions of momentum integration are invariant under boosts, we start by assigning 
the collinear subgraphs an order of magnitude ^dimension ^ contributes unity 

to the power of Q. Note that we define the collinear factors to include the integrals 
over the momenta of the loops that connect the collinear subgraphs and the hard 
subgraph. 



4 The factor 3 in Eq. (17) is the number of colors and the factor 1/2 corresponds to the spin of 
the quarks as the factor 1 in front of N g corresponds to the spin of the gluon. 
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In the next step, we have to take into account that the collinear subgraphs are 
coupled to the hard subgraph by Dirac and Lorentz indices. The effect of boosting 
a Dirac spinor from rest to the energy Q is to make its largest component of order 
(Q/m) 1 / 2 bigger than the rest frame value and the effect on a Lorentz vector is to 
give a factor of (Q/m) 1 . The exponents are just the spins of the corresponding 
fields. Multiplying by these powers gives 

Q 2 - N ". (4.9) 

This agrees with Eq. ( |4.6| ) in the case that all external lines of the hard graph are 
quarks but is a factor Q Ng larger if the external particles are gluons. 

The well-known problem of gluons with scalar polarization (see, for example, 
PH . |5"2| |) will be dealt with later on. Suffice it to say here that gluons with such a 



polarization can be factorized into the parton distributions by using gauge- invariance 
arguments. 

For the moment we just need to define the concepts of scalar and transverse 
polarization in the sense that we will use and show how this affects the power 
counting. 

Consider the attachment of one gluon, of momentum /c M , from the collinear-to- 
A subgraph to the hard subgraph. One has a factor A^g^H", where A^ and H v 
denote the collinear-to-A and hard subgraph respectively, and is the numerator 
of the gluon propagator in the Feynman gauge. On can now decompose this factor 
into components: 

A x H = A + H- + A-H + - A ± x H ± . (4.10) 
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One observes now that after the boost from the proton rest frame, the largest 
component of A^ 1 is the + component. The largest term is therefore A + H~ and this 
is the term which gives the power in Eq. (O). The other two terms are suppressed 



by one or two powers of Q. 

One can define now the following decomposition: 

A» = k^+(A»-k^y (4.11) 

The first term will be called the scalar component of the gluon and gives a polar- 
ization vector proportional to the momentum of the gluon. The second term gives 
the transverse part of the gluon with no + component, it, therefore, gives a con- 
tribution to A x H which is one power of Q smaller than the contribution of the 
scalar component. The factor in the scalar term multiplies the hard subgraph 
and gives a quantity that can be simplified by the use of Ward identities. 

The above decomposition will now be applied to every gluon joining the sub- 
graphs A and H. The contribution of our region to the amplitude is now a sum of 
terms in which each of these gluons is either scalar or transverse. Each term has a 
power 

Q 2-N q -N aQ N 3 = Q2-N q -N T . 

where N s is the number of scalar gluons and = N g — N s is the number of 
transverse gluons that enter the hard scattering. This is exactly the power in Eq. 
without a soft subgraph. 
In Taylor expanding the hard subgraph, we will need to slightly modify the above 
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decomposition in Eq. (|4.11|) and also apply an exactly analogous argument to the 
couplings of soft gluons to a collinear graph. 

One will also need to pick out the largest component of the Dirac structure of 
the collinear subgraphs, here however, this is not necessary to be made explicit since 
we do not have a cancelation of the highest power. Note that the projection of the 
largest Dirac component is directly reflected in the 7 + factor in the definition of the 
quark distribution. 

In the derivation of Eq. (|4.12| ) we assumed x not to be small. In the case of large 
x then, we have to boost parts of the collinear-to-A subgraph to get + momenta of 
order p + instead of xp + . This leads to groups of lines with very different rapidities. 
In Feynman graphs, the effect is just to give a factor 1/x times logarithms, however 



only if the exchanged lines between the regions of different rapidities are gluons |63 
In the case of exchanged quarks, there will be a suppression factor of x. All of this 
does not influence the power of Q. 

If one now adds in a soft subgraph S, one has the problem of choosing the 
appropriate scaling of the momenta. There are several possibilities and the literature 
is not entirely clear about the best scaling. In our case we will choose to have all 
soft momentum components of 0(m 2 /Q). This has the advantage of not sending 
any collinear-to-A lines far off-shell, the trade off is that we introduce regions where 
the momenta are unphysically soft in a confining theory and the power counting is 
mass sensitive. 

With this scaling the basic power in the Breit frame is m 2 /Q to a power which 
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is the dimension of the soft subgraph, including the integration over the soft loop 
momenta that circulate between S and the rest of the graph with the assumption of 
negligible masses. The power is simply 

-N gS -lN qs (4.13) 

where N g s and N q s are the number of external gluons and quarks of the soft sub- 
graph. These external lines go either into the hard subgraph in which case the 
dimension of the hard subgraph is reduced by 3/2 for each extra quark and 1 for 
each extra gluon or into the collinear-to-A subgraph which does not affect the power 
of Q. Now we have to take into account that we are dealing with vectors and spinors 
connecting S to A, this means that we have to use the same factors as before meaning 
we gain a factor Q 1 ^ 2 for each quark and Q 1 for each gluon. 

Putting all the factors together we get for the soft to hard contribution — N g — 
|iVq — N g — |iVq = —2N g — 3N q and for the soft to collinear contribution — N g — 
fiVg + N g + \N q = -N q . This together with our earlier result (Eq. ( p~T2D ) yields 
Eq. O . 

4.3.3. Leading Regions 

The leading regions for the amplitude are those with the largest exponent p(7r) in 
Eq. (|4.6|) . It is easy to see that these correspond to the reduced graphs in Fig. 
[4.3| , independently of whether the out-going photon is real or far off-shell. The 
corresponding power is Q°. These reduced graphs have direct photon couplings 
to the hard subgraph, they have exactly two parton lines connecting the collinear 
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a) b) 

Figure 4.3: Those reduced graphs that contribute to the leading regions in DVCS. 

subgraph A to the hard subgraph H, and they have no soft lines connecting to H. 
The two kinds of graph differ only by the absence or presence of a soft subgraph 
that connects to A alone. 

Among the other reduced graphs, which are non-leading for our process, are 
those of the type in Fig. [4. l|b, which are leading in the case of diffractive meson 



production, where the leading region gives Q^ 1 . 

In the case of a photon that is off-shell by order Q 2 , the amplitude for production 
of the photon behaves like Q°, the same as for a real photon. However, the physically 
observed process includes the decay of the time-like photon (to a yU + /i~ pair, for 
example), which results in a power suppression of the observed cross section by 
1/Q 2 compared with the cross section for making real photons. 
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Figure 4.4: Soft gluon loop attaching to collinear line. 



4.3.4. Proof of absence of a soft part in leading regions 



As mentioned in Sec. |4.3.3| , there might, in principle, be a soft part S in the leading 
reduced graph connected solely to the A graph by gluons, as shown in Fig. [0]b. 
Note that by Eq. (|4.6| ), quarks connecting S to A would lead to a power suppression. 
We will now show that this soft part S is indeed absent, and so we only need to 



consider regions of the form of Fig. |4.3|a. 



We will first examine a simple one loop example, Fig. |4.4j. The external quark is 



part of the A subgraph in Fig. [4.3| b, and the gluon is soft. So we parameterize the 
momenta by: 



k, 



{x x p + , k A ,k At± ), 



(4.14) 



where the k\ is 0(Q) and all the other components are of 0(m) or smaller. 

If we omit irrelevant factors in the numerator, the loop integral takes the follow- 
ing form: 

d A k 



soft k (k 2 + ie) [(kA — k) 2 — m? + it] 
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(£k 

soft k (2k+k- -k 2 ± + ie) 2(xip + - k+){k1 - k~) - (k A ,± - k±) 2 - m 2 + ie 

r dk + 

- / f ± t- 4 - 15 

Aoft k [2k+k- -k\ + ie) 2 Xl p+(k A - k~) - (k A± - k ± ) 2 - m 2 + ie 

As before, ~ stands for "equality up to power suppressed terms". As one can see, 
there is no fc + -pole in the second part of the denominator and we can freely deform 
the contour in k + to avoid the pole in the soft gluon propagator. This takes us out 
of the soft region for k. 



In the general situation, Fig. |4.3|b, we can use a version of the arguments in Ref. 



21, |65| to show that the soft momenta kf can be rerouted in such a way as to 



exhibit a lack of a pinch singularity. The essential idea is that one can find a path 
backwards or forward from one external line of S to another external line of 5*. The 
loop is completed along lines of A, all of which have much larger + momenta than 
what is typical of soft momenta, and hence there is no pinch. 

4.3.5. Subtractions 

The subtractions^] necessary to avoid double counting in the amplitude are exactly 



similar to the ones in Sec. VI of Ref. [24], since the distributional arguments to 



construct the subtraction terms on a pinch-singular surface 7r presented there are 
very general in nature and are not limited to the case of diffractive vector meson 
production that was considered in |fj| . 



Let us review the construction. For each graph T, there may be several different 
pinch singular surfaces 7r in the space of loop-momentum contributing to the leading 



5 We give examples of hands-on applications of the distributional ideas, as employed in this 
section, in the Appendix A. 
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power. We write therefore the graph as a sum of those contributions 

Asyr = EJV, (4.16) 

where 'Asy' denotes the asymptotic behavior of the graph. 

The term is obtained by Taylor expanding the hard and collinear subgraphs 
in powers of the small variables denoted by T n . Since one can possibly have several 
regions contributing to a given graph, one must make subtractions to avoid double 
counting. This operation will be denoted by R. Therefore, we can write 

Asyr = E^IV (4.17) 

This is exactly analogous to the Bogoliubov i?-operation for renormalization. 
The easiest way to formulate the above procedure is due to Tkachov et al. [p4}| . In 



this method the integrand of each graph V is viewed as a distribution, in other words 

r<g>/ = J dkT(k,p)f(k), (4.18) 

where k denotes the collection of loop momenta, p the external momenta and f(k) 
is a test function. Putting the test function to 1 will give the contribution of F to 
the amplitude. 

The advantages of the methods employed in Ref. [64"] are the control they give 



in treating different regions of momentum space separately without having sharp 
boundaries between different regions. Note that this is particularly convenient in our 
case where one has to deform contours of integration away from pinch singularities. 
If we were to use sharp boundaries these deformations would be impossible. 
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Using this language, the contribution IV to AsyT from the neighborhood of a 
pinch-singular surface tc is localized on the surface, in other words is proportional 
to a Vfunction that restricts the integration to the surface. One then defines a 
hierarchy of regions through set-theoretic inclusion: tti > 7r 2 means that the pinch- 
singular surface m contains 7r 2 . One than constructs any 7r on the assumption that 
the terms for all bigger regions have already been constructed. Thus one has a 
recursive construction starting from the largest region. 

Assume that one has constructed the terms IV for all regions bigger than ir and 
then decompose AsyT as 

AsyT = XV >7r IV + IV + other terms. (4.19) 

The "other terms" correspond to three classes of surfaces. The first class are those 
surfaces that are smaller than tt, the second class consists of those that intersect tt 
in a subset of tt, i.e. , they have a lower dimension and the third class are those that 
do not intersect tt at all. 

The hypothesis for our induction is that the sum of IV over tt' > tt gives a 
good approximation to the original T except in the neighborhoods of the smaller 
surfaces for which IV has not yet been constructed. The integrals defining the IV 's 
cover the whole space of integration variables, but are only required to give a good 
approximation if one excludes neighborhoods of smaller surfaces. This means that 
we require the test function in Eq. ( [4. 18 ) to be zero on these smaller surfaces. When 



IV is combined with the IV for larger surfaces it must give a good approximation to 
r on a neighborhood of tt. Note that a good approximation is not necessary on the 
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surfaces where the test function is zero; there is no need for IV to be constructed 
for the smaller surfaces since they will give zero with such a test function. This is 
enough to prove the inductive hypothesis for the next recursive step. 

Since IV is located on n, it is necessary to consider only a neighborhood of tt. 
This combined with the statement in the previous paragraph ensures that there is 
no need for the unconstructed "other terms" in order to construct IV 

One can therefore define 



where T w stands for a Taylor expansion in powers of the small variables on tt. The 
first term is the Taylor expansion of the original graph and the remaining terms can 
be thought of as subtractions that prevent double counting of contributions to the 
integral over a neighborhood of n. 

This results in a sum over IV and the terms for the larger regions 



giving the correct contribution to the asymptotics of T that originates from a neigh- 
borhood of 7r and of all larger regions. Neighborhoods of smaller regions are, of 
course, excluded. 

In general IV gives a divergence when integrated with a test function over a 
neighborhood of any of these smaller regions. Therefore, it is only defined when 
integrated with a test function which is on these smaller regions. We now extend 
this distribution to a distribution defined on all test functions by adding infra- 




(4.20) 



IV + xv >7r i?iv, 



(4.21) 
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red counterterms to cancel the arising divergences from the smaller regions. This 
operation is analogous to the well known + distribution employed earlier in this 
thesis. The result will be RT n with the counterterms all being local in momentum 
space. Choosing a definition of IV on the smaller surfaces is perfectly satisfactory, 
since we have not yet considered how to approximate T in regions smaller than ir. 
We only require .R1V to be finite and the counterterms to be localized on smaller 
surfaces than 7r in order not to affect the good approximation we have for n and 
larger surfaces. Continuing with the recursion, one eventually obtains appropriate 
approximations for these smaller surfaces. 

Note that the subtraction terms, as defined in Eq. ( |4.20| ), ensure that changes 
in the choice of counterterms localized on any particular surface 7r are cancelled by 
corresponding changes in the subtraction terms in the definition of IV Therefore, 
the overall result for the asymptotic expansion of V is independent of these choices. 

The above statements lead to the following asymptotic form of the amplitude 

Asy T = S r Asy T = A x H. (4.22) 

where V stands for a possible graph for the amplitude T. 
4.3.6. Taylor expansion 

We now obtain the leading term in the hard subgraph, when it is expanded in powers 
of the small momenta. The arguments used are exactly analogous to the ones used 



in Sec. VII of Ref. |24| except that we do not have to deal with a B subgraph as 
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was the case in p4 |. So we have: 

ix5 - J dk\H (q,q', (fc+0, ± ), (A+ -4, 0,0. 

/" dkJd 2 k A:± A(k A , A-k A ), (4.23) 

where fc^ is the loop momentum joining the A and H subgraphs, and again ~ means 
equality up to power-suppressed corrections. Eq. (|4.23|) has already a factorized 
form. However we still have to deal with the extra scalar gluons that may be 
exchanged between the subgraphs A and H; this will be done in the next subsection. 
Eq. ( |4.23|) can be written in the following way: 

AxH-ZiJ dk + C i (q,q',k + )O i (p,p',k + ), (4.24) 

where the Cj are the short distance coefficient functions and the Oi are the matrix 
elements of renormalized light-cone operators. 

4.3.7. Gauge Invariance 

In order to identify the Oi with the parton distributions as defined in |24| (for 
example), it is necessary to show that all gluons with scalar polarization attaching 



to the hard graph can be combined into a path-ordered exponential. Fig. 4.5 shows 



the example of one scalar gluon. We can follow the arguments of Sec. VII. D of Ref. 



24 1 step by step since those rely on very general results obtained by Collins | 66fl . In 
this way we obtain exactly the same parton distributions as in |24[, namely: 

dy 



U/v = I ^e-^ + ^(p|T^(0,y-,0 ± ) 7 + ^(0)b'), 

J-oo 47T 

00 dy~ 1 
-oo 2ir Xix 2 p~ { 



f g /p = - r T--^T e- iX2P+y ~{p\TGt(0,y-,0±)VG» + (0)\p'). (4.25) 

J— oo Z7T XiXoZ? 
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a) b) 

Figure 4.5: a) A scalar gluon attaching the collinear subgraph to the hard subgraph 
H in the unfactorized form, b) Factorized form after application of gauge invariance 
and Ward-identities. The double line represents the eikonal line to which the scalar 
gluon attaches. 

Here, V represents a path-ordered exponential of the gluon field that makes the 
operators gauge invariant. The variable X2 is the same as in Eq. ((4.3|). The evolution 



equations are the same as in [5^] and ]24] 



4.3.8. Partons with k + = 0: breakpoints and endpoints 



In the factorization theorem Eq. (|4.3| ), the integral over the fractional momenta in- 
cludes the points x% = and X2 = 0. At these points, the hard scattering coefficient 
for DVCS has a pole, and so we appear to get a logarithmic contribution to the 
cross section from a region in which one of the lines joining the parton density to 
the hard scattering subgraphs is soft instead of collinear. This apparently contra- 
dicts our power-counting result that such a region gives a non-leading power. This 
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phenomenon was investigated by Radyushkin In this section, we will give a 

general demonstration that the region in question does not give a problem. 

First, let us observe that the region of integration over x\ in the factorization 
formula Eq. ( f4.3| ) is— 1 + x < X\ < 1. This is proved by the methods of light- 
front perturbation theory, by requiring that the intermediate states in Fig. [4.3| be 



physically allowed. See Ref. pi ffTJ for detailed derivations and discussions. The 



points x\ = and X2 = at which the potential problem arises are what we will 
call "breakpoints" , since they occur in the middle of the range of integration where 
one of the two lines changes direction. 



We continue by examining a particular case, illustrated in Fig. |4.6| , and showing 
how the argument generalizes. To simplify the example, let us restrict our attention 
to regions where the subgraph A and the lower three lines (p, p' and k — p') have 
their momenta collinear to the proton. We will also require the two quark lines, 
k + p — p' and k, on the sides of the ladder to have their momenta either collinear 
to the proton or soft. 

We will also only need the case of the production of a real photon, q' 2 = 0, since 
this is where the problem arises. 

The top loop of the graph has the form, omitting the ie factors in the propagators 
for convenience : 

U [ d A k Numerator factors ,^ ^ 

J {k 2 — m?) \{k — p') 2 — m 2 ] \{k + p — p') 2 — m 2 ] \{k + q') 2 — m 2 ] 

When both k and k + p — p' are collinear to A, the top line is off-shell by 0(Q 2 ), 
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Figure 4.6: Particular example of potentially problematic diagram, 
and it is correct to use the collinear approximation 

\ ► I ; (4.27) 

(k + q') 2 — m 2 + ie X2Q 2 /x + ie 

where £2 = k + /(xp + ). A corresponding replacement is also to be made in the 
numerator in U. The right-hand side of Eq. (^4.27j) exhibits the afore-mentioned 
pole at X2 = 0. The result of applying the collinear approximation is to give the 
appropriate contribution to the factorization formula Eq. ( f4.3| ). 

The collinear approximation becomes invalid when k becomes soft, i.e., when 
X2 — > 0. We must now demonstrate two facts. The first is that, when k is in a 
neighborhood of the soft region, the collinear approximation is valid after integration 
over k. The second fact is that the use of the collinear approximation does not give 
an important contribution from some other region of k that was not permitted in 
the original graph. 
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We now examine the integral U in the neighborhood of the soft region for k. It 
has the following form, again omitting the ie factors in the propagators : 

Usoitk ~ Lit k ^ [2k+k- -k\- m 2 ] [2p' + (p'- - k-) - (p> - k)l - m 2 ] (4 ' 28) 

1 

[2(p+ - p' + )(p- - p'~ + k~) -(p-p' + k)\ - m 2 } -k 2 - m 2 ) ' 

where we have neglected k + in the collinear-to-A lines, k~ in the collinear-to-S line 

and x ^ x\. We have also ignored the numerator factors, which are an irrelevant 

complication for our purposes. 



According to the power-counting results of Sec. 4.3.2, which are obtained from 



p5 |, the soft region for k gives a power-suppressed contribution. This estimate as- 



sumes that all components of k are comparable (in the Breit frame), and is obtained 
as follows. Let the magnitude of the components of k be m. Then the order of mag- 
nitude of the soft part of U is a product of factors l/(m 5 Q 3 ) from the denominators, 
m 4 from the phase space, and Q 2 m 2 from the numerator, for an overall power m/Q. 
This result can be obtained by writing down the largest components in the trace 
and propagators of Fig. 16 and Eq. ( 4.28Q . Moreover in this region it is correct 



to replace the fourth propagator in Eq. (|4.28| ) by its collinear approximation Eq. 
( |4.27| ), so that we do not lose the factorization theorem. 

However, if the components of k are asymmetric this estimate no longer holds. 
In particular if the longitudinal components of k are of order k^ ~ m 2 /Q while 
the transverse components remain of order m, then we get contributions of order 
1/m 8 from the denominators, m 6 /Q 2 from the phase space, and Q 2 m 2 from the 
numerators, for a total of m°Q°. This shows that the contribution from this region 
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is unsuppressed for large Q. 



At this point we must appeal to the contour deformation arguments of Ref. |65 
It is only when the integration over k is pinched in the region in question that it 
needs to be taken into account. In the dangerous region we have k + k~ -C kj_, so 



that the only k + dependence in Eq. ( f4.28| ) is the pole in the fourth denominator. 



We can therefore deform k + into the complex plain a long way out of the region we 
are considering, indeed all the way to the collinear-to-A region. Then the collinear 
approximation is valid so that we can replace the graph by its contribution to the 
factorization formula. 

This contour deformation argument is completely general, as explained in Sec. 
HIE of Ref. H^J. Whenever we have a soft momentum with k + k~ <C kj_, the contour 
of k + can be deformed away from poles in the jet subgraph associated with with the 
produced photon. Since all the relevant singularities are in the final state, they are 
all on the same side of the real axis. 

Now that we have established in more detail that the only leading regions are 
those symbolized in Fig. [|7|(a), we can apply the collinear approximation as de- 
scribed earlier, and hence we obtain the factorization theorem. 

But we still see the following problem. In the factorization theorem, Eq. ( |4.3[ ), 
the parton densities are non-analytic at the breakpoints x\ = and 22 = 0, whereas 
the coefficient function has a pole at each of these points. Again consider the 



collinear approximation to Fig. |4.6| in the region we were considering. The parton 



density is non-analytic when X2 = 0, while the coefficient function has a pole there, 
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as is seen from the right- hand- side of Eq. (|4.27 ). So we cannot literally apply the 
contour deformation argument. 

What we will show is that the parton density is continuous at the breakpoint, 
so that it can be written as the sum of a function that is analytic at x<i = and a 
function that has a zero at x<i = 0. The only potential leading twist contribution 
near the breakpoint is associated with the non-zero analytic term to which the 
contour deformation argument applies. 

To prove this property of the parton density at a breakpoint, consider a general 



graph for the parton density, as shown in Fig. [4.7| . We have found it convenient to 
change the labeling of the momentum compared with the previous figure. As always, 
the k~ and components of k have been short circuited and are integrated over. 
The /c-line gives a pole at k~ = (k\ + m 2 — ie)/2k + = (k\ + m 2 — ie) /2x\P + , while 
the k + q' — g-line gives a pole at k~ = [(k + q — q')\ + m 2 — ie]/2(k + — £P + ) = 
[(k + q — q')\ + m 2 — ie]/2x2-P + . Here, £ is the fractional longitudinal momentum 
transfer 1 — P' + / P + . In addition there are poles from the collinear-to-A lines in the 
blob. For example if the blob consists of a single line, we have a pole at k~ = P~ — 
{k 2 + m 2 ± -te)/2(l-x 1 )P + or at Ar = -P'~ + [(k + P') 2 + m 2 ± -te)/2(l-£ + x 1 )P + . 

As we vary xi, the k~ contour can generally be deformed to avoid the poles, so 
that we have analytic dependence on x\. The possible exceptions occur when the 
k~ contour is pinched for finite k~ or when a singularity coincides with the endpoint 
of the integration at k~ = oo. A pinch never occurs; in the general case this is a 
consequence of the Landau rules. But endpoint singularities occur, and these are 
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precisely at the breakpoints. 

For example if k + — > 0, then k + can approach from above and below. The pole 
giving us trouble stems from the /c-line, all other propagators are unproblematic in 
this case, since their poles are at finite k~. The pole in k~ approaches +00 — ie as 
one approaches + and — 00 -Me as one approaches 0~. This means that the k~ pole 
crosses the real axis at infinity. Hence the parton distribution is non-analytic there. 
Since the singularity is at \k~\ = 00 the other propagators have large denominators, 
and hence we get a zero for the non- analytic part of integral at the breakpoint. Thus 
the parton density is continuous at the breakpoint, as claimed. This result enables 
the factorization formula to be valid in the neighborhoods of the breakpoints. Since 
the other poles in the k~ integral are on opposite sides of the real axis, the parton 
distribution is non- zero at the breakpoints. 

Effectively the crossover of the pole occurs when A; is in a collinear-to-i? region, 
which we know is power suppressed. This indicates that the argument we have just 
given generalizes to all graphs. 

We also remark on the behavior at the endpoints. Let us look at the case 
k + — > p + . We find that another of the poles pole runs off to —00 this time and 
crosses the real axis there. But now all the other poles are on a single side of the 
real axis, so that the sole contribution to the parton density comes from the pole at 
infinity, and hence there is a zero of the parton density at the endpoint. 
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Figure 4.7: Parton distribution amplitude. 
4.3.9. Completion of Proof 

Using the definitions of the parton distributions and the hard scattering coefficients 



a real photon which directly goes into the final state and for the production of a 
time-like photon that decays into a lepton pair. 

4.4. Conclusion 

We have proved the factorization theorem for deeply virtual Compton scattering up 
to power suppressed terms to all orders in perturbation theory. The form of the 
theorem is independent of the virtuality of the produced photon. 



we finally obtain Eq. ( [4.3| ). Note that the theorem is valid for the production of 
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CHAPTER 5 

Diffractive Exclusive Photon Production in DIS at HERA 



5.1. Introduction 

Recent data from HERA has spurred great interest in exclusive or diffractive direct 
production of photons in e — p scattering (deeply virtual Compton scattering or 
DVCS) as another source to obtain more information about the gluon distribution 
inside the proton for nonforward scattering. Therefore, we will start a phenomeno- 



logical investigation of this process in this chapter based on Ref. [bo after having 
proved the validity of a factorization theorem for this process in the last chapter. 

Exclusive, diffractive virtual Compton processes at large Q 2 , first investigated 
in [^7[ , offer a new and comparatively "clean"[] way of obtaining information about 
gluons inside the proton in a nonforward kinematic situation. We are interested 
in the production of a real photon compared to the inclusive DIS cross section. 
The exclusive process is nonforward in its nature, since the photon initiating the 
process is virtual (q 2 < 0) and the final state photon is real, forcing a small but 
finite momentum transfer to the target proton, i.e. , forcing a nonforward kinematic 
situation as we would like. 



1 Clean in the sense that the wave function of a spatially small size configuration within a real 
photon is better known as compared to the wave functions of vector mesons thereby removing a 
big theoretical uncertainty in the determination of the gluon distribution. 
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The chapter is organized in the following way: In Sec. |5.2| we estimate the 
amplitude in the normalization point Q$ ~ 2GeV 2 using the aligned jet model 
approximation and conclude that for such Q 2 the nondiagonal amplitude is larger 



than the diagonal one by a factor of ~ 2. In Sec. |5.3| we calculate the imaginary 
part of the amplitude for 7* +p — > 7 + p in the leading order of the running coupling 
constant a s and compare it to the imaginary part of the amplitude in DIS in the 



same order. In Sec. |5]^ we argue that at sufficiently small x the t-dependence of the 
cross section should reflect the interplay of hard and soft physics typical of diffractive 
phenomena in DIS. In other words, at fixed x and increasing Q 2 , hard physics should 
tend to occupy the dominant part of the space of rapidities. In contrast to this, at 
fixed Q 2 and decreasing x, hard physics should occupy a finite range of rapidities 
which increases with Q 2 as ~ In -0L anc i with f3 ~ 0.1 — 0.2 in the HERA energy 
range due to the QCD evolution, and that soft QCD physics occupies the rest of 



phase space. In Sec. and Sec. ^]6] we give the total cross section of exclusive 
photon production and give numerical estimates of the DVCS production rate at 
HERA and find that such measurements are feasible for the current generation 
of experiments. We also show the feasibility of directly measuring the real part 
of the DVCS amplitude and hence, at least, the shape of the nondiagonal parton 
distributions through a large azimuthal angle asymmetry in ep scattering for HERA 
kinematics. Sec. |5.7| finally contains concluding remarks. 
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5.2. The amplitude for diffractive virtual Compton scattering at interme- 
diate Q 2 



Similar to the case of deep inelastic scattering, in real photon production it is possi- 
ble to calculate within perturbative QCD the Q 2 evolution of the amplitude but not 
its value at the normalization point at Qq ~ few Ge V 2 where it is given by nonpertur- 
bative effects. Therefore, we start by discussing expectations for this region. It was 
demonstrated in that the aligned jet model |J7| coupled with the idea of color 



screening provides a reasonable semiquantitative description of F 2 n(x < 10^ 2 ,Ql). 
In this model the virtual photon interacts at intermediate Q 2 and small x via tran- 
sitions to a qq pair with small transverse momenta - ko,t ((^o,t) ~ 0.15GeV 2 ) and 
average masses ~ Q 2 which thus carry asymmetric fractions of the virtual photon's 



longitudinal momentum. Due to a large transverse color separation, b ~ 2y2/3rv, 
the aligned jet model components of the photon wave function interact strongly with 
the target given by the cross section a to t{ 1 ' A J M" — N) ~ a to t{^N). Neglecting con- 
tributions of the components of the 7* wave function with smaller color separation, 
one can write 0" tot (7*iV) using the Gribov dispersion representation |69| as p*| : 



3tt J mS (Q 2 + M 2 ) 2 

where the factor M 2 in the numerator is due to the overall phase volume, 
R e+e (M 2 ) = a ^ e l e Z*^ a ^°^ . The factor 3 ^° 2 ^ is the fraction of the whole phase 
volume occupied by the aligned jet model , and the factor 1/(Q 2 + M 2 ) 2 is due to the 
propagators of the photon in the hadronic intermediate state with mass square equal 



to M 2 . Based on the logic of a local quark-hadron duality (see e.g. f7C§ and refer- 
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ences therein), we take the lower limit of integration to be Mq ~ 0.5GeV 2 < m 2 . In 
the case of real photon production, the imaginary part of the amplitude for t = is 
given by 



s w ' yt_u 3tt7m 2 (Q 2 + M 2 )M 2 

(5.2) 

with s = 2qQrriN being the flux factor. The only difference from Eq. |5.1| for <Jtot{l* + 
N) is the change of one of the propagators from 1/(Q 2 + M 2 ) to 1/M 2 - here go is 
the energy of the virtual photon in the rest frame of the target. 

Approximating R e+e (M 2 ) as a constant for the Q 2 range in question [] we find 



K ~ ImA(r + N - 7 + N) t=0 -Qi + M 2 ln {1 + Q ^ ] 



In the following analysis we will take Q\ for the perturbative QCD evolution as 
2.6 GeV 2 to avoid ambiguities. It is easy to convince oneself that for Mq 2 ~ 0.4 
0.6 GeV 2 and Q 2 w 2 - 3 GeV 2 Eq. (|J) leads to i2 w 0.5. A similar value of R 



has been obtained within the generalized vector dominance model in Ref. [71] As 
we will see below, QCD evolution leads to a strong increase of Q 2 ImA{^* + — > 
7 + N) t =o for increasing Q 2 and fixed x. However, this does not change the value 
of R appreciably. 

5.3. The amplitude for exclusive real photon production at large Q 2 . 

The process of exclusive, direct production of photons in first nontrivial order of 

a s ln %r at small x can be calculated as the sum of a hard contribution calculated 
s Qo 

2 We understand this in the sense of a local duality of the hadron spectrum and the qq loop. 
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(see Fig. |5.1| ) within the framework of QCD evolution equations [gfj and a soft con- 
tribution which we evaluated in the previous section within the aligned jet model. 
The hard contribution can be described through a two gluon exchange of the box 
diagram with the target proton. In order to calculate the imaginary part of the am- 
plitude, we need to find the hard amplitude from the box as well as the gluon-nucleon 
scattering plus the soft contribution at Qq from the aligned jet model analysis. This 
is an example of the usual subtraction scheme in factorization where one subtracts 
out the collinear region in the factorization formula to have a convolution of a par- 
ton distribution with a hard scattering coefficient and then adds it back in which 
corresponds to the aligned jet model contribution in our case. 

Let us first give a general expression for the imaginary part of the amplitude and 
then proceed to deal with the gluon-nucleon scattering, followed by the calculation 
of the box diagrams. 

First, let us discuss the hard contribution which actually dominates in the con- 
sidered process. To account for the gluon-nucleon scattering, we work with Sudakov 
variables for the gluons with momenta p\ and P2 attaching the box to the target 
and the following kinematics for the gluon-nucleon scattering: 

Pi = aq' + x\p' + p t , d 4 pi = -dadxid 2 pt, (5.4) 

where q' and p' are light-like momenta related to p, q the momenta of the target 
proton and the probing virtual photon respectively, by: 

q = q' - X p, p = p' + — -q', 

2p' ■ q' 
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Figure 5.1: Leading contribution to DVCS at small x. 

s = 2p ■ q = 2p' ■ q' — xp 2 , (5.5) 

with x being the Bjorken x and X\ the proton momentum fraction carried by the 
outgoing gluon. Note that p' is not related to P', the momentum of the outgoing 
proton. Equivalent equations to Eq. ( |5.4|) apply for pi with the only difference being 
that X\ is replaced by X2, the momentum fraction of the incoming gluon, signaling 
that in leading log there is only a difference in the longitudinal momenta but not 
in the transverse momenta. Furthermore there is a simple relationship between x% 
and X2 : A = x\ — X2 = const., following from the kinematics of the considered 

reaction]] where A is the asymmetry parameter or skewedness of the process under 

3 In the case of the imaginary part of the amplitude which we are discussing at this point, one 
has xi > A > and we can treat the soft part as a parton distribution function (the DGLAP 
regime), whereas if < x\ < A one would have the situation of a distributional amplitude as first 
discussed by Radyushkin [E3| which is governed by the Brodsky-Lepage evolution equations. 
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consideration. Therefore, one is left with just the integration over p\ since p 2 cannot 
vary independently of p%. Being exclusively interested in the small x region, one 
can safely make the following approximations: s = 2p ■ q ~ 2p' ■ q' and p' ~ p. 
Since we are working in the leading a s In Q 2 approximation, neglecting corrections 
of order a s , the main contribution comes from the region p 2 << Q 2 , hence the 
contribution to the imaginary part of the amplitude simplifies considerably. First, 
since \p\\ = \ax\s + p 2 \ << Q 2 , one has a « 1 and the polarization tensor of the 
propagator of the exchanged gluon in the light-cone gauge q'^A* 1 = becomes, see 



Ref. [31 



d,x * 44- ( 5 - 6 ) 

In other words it is enough to take the longitudinal polarizations of the exchanged 
gluons into account. 

Using Eq. (|5.6| ) one obtains the following expression for the total contribution 
of the box diagram and its permutations: 

ImA = f ^^^2ImAfPlmAf} T) dMd UCT (p 2 ), (5.7) 
J [Air) ipiP 2 

where ImA a ^ p ^ = ImA a ^ v (^f*g — > qq) is the sum of the box diagrams, ImA ^^ 
is the amplitude for the gluon-nucleon scattering, a,b are the color indices and 
the overall tensor structure has been neglected for now. The usage of the imaginary 
part of the scattering amplitude and in particular limiting ourselves to the s-channel 
contribution as the dominant part in both the forward and the nonforward case (Eq. 
(|5.71 )) is correct (see Ref. for more details) as long as we restrict ourselves to 
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the DGLAP region of small x and thus small t, where t — {jp\ — P2) 2 is the square 
of the momentum transfered to the target. The real part of the amplitude will be 
evaluated below by applying a dispersion relation over the center of mass energy s. 
Using Eq. ( |5.6| ) and the following Ward identity which is the same as in the abelian 
case since the box contains no gluons, i.e. , is color neutral: 

< (P W = 0, <j P W = 0, (5.8) 
yielding in the expression of the imaginary part of the amplitude Eq. 



ImAfPv^y = ImA a ^p tfl p tu 
4(jo • q) 2 X1X2S 2 

one can rewrite Eq. (|5.7|) as: 



(5.9) 



1mA f 1 dxi x A 2 2 2 f sdad 2 p t 2 ±ImA a £pq x q (T 

= / E{— , — ,Q ,p t ,Qo) / 2 2 Pt^a - 2 , (5.10) 

s Jx x\ x\ x\ J [Ztt] ) P1P2 s z 

where we have used < pt^Ptv >= (the average over the transverse gluon 

polarization) and defined the imaginary part of the hard scattering to be given by: 

E(^A,Q 2 ,P 2 t ,Q 2 o) = -k„ »—5 ab , (5.11) 

x\ X\ 2 M x 2 s 

where the sum over repeated indices is implied. Up to this point, we have just 
rewritten the equation for the imaginary part of the total amplitude but have not 
identified the different perturbative and non-perturbative pieces. In the case of a 
virtual photon with longitudinal polarization, this would be an easy task since the qq 
pair would only have a small space-time separation and we could follow the argument 



in Ref. |21], |24], |70[ stating that the box is entirely dominated by the hard scale Q and 
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thus can unambiguously be calculated in pQCD. However, in our case we are dealing 
with a virtual photon which is transversely polarized and thus one can have large 
transverse space separations between q and q. The resolution to this problem can 
be found in the following way: one accepts that one has a contribution from a soft, 
aligned-jet-model-type configuration and that there is no unambiguous separation 
of the amplitude in a perturbative and non-perturbative part up to a certain scale 
Qq. However, in the integration over transverse gluon momenta, one will reach a 
scale at which a clear separation into perturbative and non-perturbative part can 
be made and hence we can unambiguously calculate albeit not the imaginary part 
of the amplitude of the upper box but its In Q 2 derivative, i.e. , its kernel convoluted 
with a parton distribution. At this point then, one can include the non-perturbative 
contribution of the aligned jet model into the initial distribution of the imaginary 
part of the total amplitude and solve the differential equation in Q 2 . One obtains 



the following solution for the imaginary part |26 



ImA(x,Q 2 ,Q 2 )=ImA(x,Q 2 )+ ^_ / —P qg (—, —)g{x x , x 2 , Q' 2 ), (5.12) 

JQq Jx X\ X\ X\ 

where P qg is the evolution kernel^ and starting from Qq the gluon distribution can be 
defined from Feynman diagrams in the leading a s In Q 2 approximation by realizing 
that in Eq. Q5.10D one can replace p\ and p\ by p\ and one finds: 



"'^. ^gW = 9{X1 , X2 ,Q% (5.13) 



(2 T )V? " 



4 In Ref. ]23|, Q a similar equation was derived for the complete amplitude for larger x ~ 0.1, 
where the quark distribution dominates and one only needs the P qq kernel. Of course, at sufficiently 
small x the contribution of this term is numerically small. 
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where g is the nondiagonal parton distribution in general. Comparison of Eq. (|5.10 ) 



with the QCD-improved parton model expression for the total cross section of charm 



production given in shows that g in the case A = is the conventional, diagonal 
gluon distribution. 

Note that the parton distribution which serves as an input in Eq. ( |5.12 ) has 
to be evolved over the Q 2 -range covered by the Q' 2 integral which complicates the 
calculation. We will explain below how to deal with this issue in practical situations. 

At this point we would like to comment on equivalent definitions of nondiago- 
nal parton distributions in the literature which differ by kinematic factors (see for 
example |24|, |25|, |23|, |2|]). Eq. (|5.13| ) corresponds to the definition used in p3| [42 



since it is given on the level of Feynman diagrams. For the non-perturbative input, 
ImA(x, Ql) we will be able to use the aligned jet model analysis of Sec. |5]2| and the 
standard relation between ImA 1 *^ 1 p (x, Q 2 , t = 0) and F 2p (x, Q 2 ): 

ImA^ p ^* p (x, Q 2 , t = 0) = F2p } X : Q2) . (5.14) 

Following the discussion above, we now only need to calculate P qg to leading 
logarithmic accuracy, in order to make predictions for the imaginary part of the 
whole amplitude. Therefore, let us now consider the box diagram where the two 
horizontal quark propagators are cut, corresponding to the DGLAP region, i.e. , 
neglecting the u-channel contribution. 



The kinematics (see Fig. |5.2j) for the calculation of the cut box diagram, again 



using Sudakov variables^, is the following. The quark-loop momentum k is given 
5 Note the slight differences between the Sudakov decomposition used here and before e.g. x\ — > 
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Figure 5.2: Cut box diagram giving the kernel for the imaginary part of the DVCS 
amplitude. 

by: 

k = aq' + (3p' + p t , d 4 k = -dad(3d 2 k t , (5.15) 
where q' and p' are light-like momenta related to p, q by: 

q = q -xp p 1 =p + — - q , 

2p ■ q 

s = 2p x q = 2p ■ q' - xp\. (5.16) 

The momenta of the exchanged gluons, in light cone coordinates, are given by: 

Pi = (x 1 p + ,m 2 /Q,p t ), p 2 = (x 2 p + ,m 2 /Q,p t ). (5.17) 

The probing transverse photon and the produced photon have the following mo- 
(3 etc. 
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menta, again in light cone coordinates: 



'=<-^2^ '>- ^ -2^°'>- (5 ' 18) 
P qg is calculated in the light cone gauge yielding the following result for the most 
general caseQ: 

„ ,x A. , 9 a s „ T x(x — A) + (x\ — x) 2 

p ^'^ = *^ a i Nr • (5 - 19) 



The DIS kernel is analogous to Eq. (|5.19 ) except that A = and the kernel for real 



photon production is obtained for A = x. 

We now can proceed to calculate the total imaginary part of the amplitude from 
Eq. ( |5.12| ) where we parameterize the gluon distribution at small x as 

g(x 1 ,x 2 ,Q 2 )=A (Q 2 )xf liQ2) . (5.20) 

We neglect the x 2 dependence for the moment]] . The above parameterization is 
taken from CTEQ3L as well as the parameterization of a s in terms of Q 2 in leading 
order: 

MQ 2 ) = exp[-0.7631- 0.7241 In fi^)-1171n 2 fi^) 



+0.534 In 



/ ln(Q/A) V 
^ln(g /A)/ J 



M Q 2 ) = -0.3573 + 0.3469 In 0.3396 In 2 

lln(Q /A)y \\n(Q /A)J 



. f ln(Q/A) \ 



+ °-° 91881n ih&J' (5 ' 2i: 



6 Note that this expression is defined differently from the gluon — > quark splitting kernel as 
given in e.g. Ref. [^3], Q by a factor of 1/xi due to the fact that the additional X\ already appears 
in the convolution integral for the In Q 2 derivative. 

7 This effect will be taken into account in the actual numerical calculation - see discussion below. 
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with A, Q and a s given by: 



A = 0.177GeV, g = 1.6GeV, a s 



4tt 



91n(Q 2 /A 2 ) 



(5.22) 



where we have taken N c = 3 and N F = 3. 

The ratio R of the imaginary parts of the amplitudes^ is given by: 



R 



ImA{^f* + p — > 7* + p) 
ImAi^* + p — > 7 + p) 



(5.24) 



We give R in the x range from 10 -4 to 10 -2 and for a Q 2 of 3.5, 12 and 45 GeV 2 
since this kinematic range is relevant at HERA. One might ask, what about the 
contributions due to quarks. The answer is that the corrections are smal!0 but for 
completeness we include them here. Eq. ([5.12 ) is then augmented with a similar 
expression for the quark contribution where the kernel is now that of quark-quark 
splitting and the nondiagonal parton distribution is that of the quark: 



ImA(x,Q 2 ,Q 2 ) = ImA(x,Ql)+ f 



Q 2 dQ' 2 fdxi 
Ql Q 



x X\ 

+P qq (— ,— )q{x 1 ,x 2} Q' 2 )] 

X\ X\ 



I ^[P q9 {-A)g{^x 2l Q' 2 ) 

Jx Xi X\ X\ 

(5.25) 



where the general expression for the kernel, after a similar calculation as before, is 
found to be: 

_x x 3 A_ ( _x i_ x 2 ' 

X\ ~xf X\ \X\ ~x^- 



— = ^ 2 a—C F 



X\ X\ 



7T 



^The tensor structure which is the same in both cases, namely: 



p-q 

cancels out in the ratio! 

9 We found them to be around 10% in the ratio R. 



p-q 



(5.26) 



(5.23) 
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and the + - prescription is the one used in Ref. [41]. The quark distribution itself 
is also taken from CTEQ3LQ and given by: 

q(x u x 2 ,Q 2 ) = A (Q 2 )xf\ (5.27) 

with 

A (Q 2 ) = exp[0.1907 + 0.04205 In 0.2752 In 2 

\HQ /A)) \\n(Q /A)J 

-0 3171 In 3 f MQ/A) V 
A 1 = 0.465. (5.28) 

We chose A\ to be constant since it varies only between 0.4611 and 0.468 in the 
Q 2 range of interest, i.e. , the error we make is almost negligible since the quark 
distribution themselves are small in the x-range considered. Furthermore, according 



to our discussion in Sec. |5^, we chose the initial distribution for the imaginary part 
of the DVCS amplitude to be twice that of the initial distribution for the imaginary 
part of the DIS amplitude. In the evolved QCD part, the nonforward kinematics are 
taken into account in the kernels of the QCD evolution equation, also the different 
Q 2 evolution of nondiagonal as compared to diagonal distribution has been taken 
into account as explained below. 

As the calculation with MATHEMATICA showed, the amplitude of the produc- 
tion of real photons is larger than the DIS amplitude over the whole range of small 
x and R turns out to be 0.551, 0.573 and 0.57 for x = 10~ 4 , 0.541, 0.562 and 0.557 



10 We used the u-quark parametrization for all light quarks for simplicity, which is surely un- 
problematic at small x. 
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for x = 10~ 3 and 0.518, 0.519 and 0.505 for x = 10~ 2 in the given Q 2 range. It has 
to be pointed out that for a given Q 2 , the ratio is basically constant. Of course, the 
ratio R will approach 1/2 as Q 2 is decreased to the nonperturbative scale since this 
is our aligned jet model estimate. The reason for the deviation from R— 1/2 is due 
to the difference in the evolution kernels. 

It is worth noting that in the kinematics we discuss, the ratio is still rather 
sensitive to the nonperturbative boundary condition. For example, assuming the 
same boundary conditions for DVCS and DIS, would result in a reduction of R of 
about 20(10)% at Q 2 ~ 12(40) GeV 2 and x ~ 10~ 3 

In Eq. ( |5T25|) the median point of the integral, as found by a MATHEMATICA 
program, corresponds to x±/2 ~ x 2 ~ x. This is due to the mass of the qq pair in the 
quark loop being oc Q 2 . For such a Xi/x 2 the ratio of nondiagonal to diagonal gluon 
density depends weakly on x 2 . Therefore, with an accuracy of a few percent, we can 
approximate this ratio by its value at X\jx% = 2. Therefore, in the calculation of R, 
we used Eqs. ( |5.20| ) and fl5.27| ) for both the diagonal and nondiagonal case but then 



multiplied the real photon result of the amplitude by a function f(x, Q 2 ) for each x 
and Q 2 to take into account the different evolution of the nondiagonal distribution 
as compared to the diagonal one, 

t At r>2 rfl\ t A( rt\ , [ Q * dQ' 2 f(x,Q' 2 ) ^ dx x ( x A 2 
ImA(x, Q ,Q ) = ImA(x, Q )+ 7^ / P qg {— , — )9{xi, x 2 , Q )■ 

JQl Q' z Jx X\ X\ Xi 

(5.29) 

The function was determined by using our modified version of the CTEQ-package 
and, starting from the same initial distribution, evolving the diagonal and nondiag- 
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onal distribution to a certain Q 2 and comparing the two distributions at the value 
X2 = x\/2 = x for different x and then interpolating for the different ratios of the 
distribution in Q 2 for given x. For this median point the difference between the 
diagonal and nondiagonal gluon distribution is between 8 — 25% depending on the 
x and Q 2 involved and 10 — 55% for the quarks (see the figures in Ref. |3(J for 
more details). 

As far as the complete amplitude at small x is concerned, we can reconstruct 
the real part via dispersion relations |32| |33|, which to a very good approximation 



gives: 

Re A 7r d In(xImA) 

"=7^1=2^1^' (5 ' 30) 

x 

Meaning that since 1mA can be fitted as x _1_<5 , r\ ~ |5 is independent of x to a good 
precision. Therefore, our claims for the imaginary part of the amplitude also hold 
for the whole amplitude at small x. This is due to the fact that within the dispersion 
representation of the amplitude over x, the contribution of the subtraction constant 
becomes negligible at sufficiently small x. 

One also has to note that there is a potential pitfall since the QED 
bremsstrahlung - the Bethe-Heitler process - , where the electron interacts with 
a proton via a soft Coulomb photon exchange and the real photon is radiated off 
the electron, can be a considerable background. As was shown by Ji pS ], the Bethe- 



Heitler process will give a strong background at small t and medium Q 2 and x > 0.1. 
We will discuss this subject in more detail later on. 
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5.4. The i-slope of the 7*iV — > 7iV cross section 

The slope of the differential cross section of the virtual Compton scattering 
daJ dt ^ ex P(Bt) is determined by three effects: (i) the average transverse size 
of the qq component of the 7* and 7 wave functions involved in the transition, (ii) 
the pomeron-nucleon form factor at the nucleon vertex, and (iii) Gribov diffusion 
in the soft part of the ladder. This leads to several qualitative phenomena. At the 
normalization point, qq configurations of an average transverse size, comparable to 
that of the p-meson, give the dominant contribution to the scattering amplitude, 
leading to a slope similar to that of the processes 7 + p — > p,u + p. The contribution 
of the higher mass qq components is known to result in an enhancement of the dif- 
ferential cross section of Compton scattering at t = by a factor pa 2 as compared to 
the prediction of the vector meson dominance model with p,u, <f), J/ip intermediate 
states, see e.g. |jT5| . Since the diffraction of a photon to masses Mx > l.ZGeV 
has a smaller t slope than for transitions to p and u, one could expect that the 
high mass contribution would lead to a t-slope of the Compton cross section being 
somewhat smaller than for the production of p, cu-mesons. However direct experi- 



mental comparison |73| of the slopes of the Compton scattering and the a;-meson 
photo-production at (Ej nc ) pa lOOGeV finds these slopes to be the same within the 
experimental errors. Using these data, we can estimate the slope of the amplitude 
for diffractive photon production in DIS at HERA energies but at moderate Q-i.e. 
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in the normalization point as 

B(s,Ql) = Bcom P .Scatt.(so) + 2c/ln(— ), (5.31) 

So 

where a' = 0.25Ge1/- 2 , s = 200GW 2 , and B Comp .scatt.{ s o) = 6.9 ± 0.3GeV- 2 Q 
Hence for HERA energies B(W = 200GeV r , Q 2 ) ~ lOGel/" 2 . 

In another limit of large Q 2 and large enough x, say a; ~ 10 -2 , the dominant qq 
configurations have a small transverse size and the upper vertex does not contribute 
to the slope. Furthermore, the perturbative contribution occupies most of the ra- 
pidity interval and leaves no phase space for the soft Gribov diffusion. In this case, 
the slope is given by the square of the two-gluon form factor of the nucleon which 
corresponds to B = B ggN m 4 5GeV~ 2 [ pTf . 

An interesting situation emerges in the limit of large but fixed Q 2 when the 
energy starts to increase. In this case, the perturbative part of the ladder has the 
length ~ ln(^^). Here k = x/xq, where xq is the fraction x of the parent parton 
at a soft scale. For HERA kinematics k ~ 0.1 - 0.3 for Q 2 ~ 10 - 20GeV 2 and 
decreasing with increasing s. Thus at high Q 2 one has an approximate factorization 
for diffraction in the case of high masses (M 2 > lOOGeV^ 2 , M 2 ^> Q 2 ) in the 
scattering of real and virtual photons observed at HERA [74], namely 



1 da(-?N -> XN)(W, M x ) _ 1 da(YN XN)(W, M X ) ^ ^ 



a^N) dtdM 2 VtotiTN) dtdM 2 

The observed slope for these processes is B ~ 7GeV~ 2 which is consistent with the 

presence of a cone shrinkage at the rate ~ 2a' ln(W 2 /M 2 ) as compared to the data 

11 Note that the data J73| can be equally well described by the fit da/dt oc exp(Bt) with B = 
6.9 ± 0.3GeV~ 2 and by the da/dt oc exp(8.9i + 2.2t 2 ) fit. 
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at lower energies where smaller values of W 2 /M 2 were probed. Similarly we can 
expect that for virtual Compton scattering at large Q 2 , the slope will increase with 
decrease of x, at very small x, approximately as 

B(W 2 , Q 2 ) QW = B ggN + 2a'(HW 2 K/Q 2 ) - ln(W 2 /m 2 p ))9(W 2 n/Q 2 - W 2 /m 2 p ), 

(5.33) 

where Wq = 200GeV 2 . We take into account here that B 99 n was determined ex- 
perimentally from the processes at W 2 ~ W 2 . 

5.5. The rate of exclusive photon production at HERA 

To check the feasibility of measuring a DVCS signal against the DIS background, 
we will be interested in the fractional number of DIS events to diffractive exclu- 
sive photoproduction events at HERA in DIS which will tell us whether it will be 
statistically feasible to search for DVCS events among DIS events. We define this 
fractional number of events as: 

#7 = TV^ T + *=o x ~5 a tot{l PJ 5.34 

°tot\l P) at B 

Using 

^(7* + P - 7 + P) = + V 2 )e B \ (5.35) 

which can be derived from applying the optical theorem and using R, the ratio of 
the imaginary parts of the amplitudes given by Eq. ( [5.24 ), rj = Re A/ Im A as given 



by Eq. ( |5.30|) and where t = —-jz^. Pt — ~p\ with t min = — ~ 0, one can 



now rewrite Eq. ( |5.34j) . A complete expression for DVCS will be given in the next 
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section. Note that only da ^ + ^ 1+p ^ \t=o is calculable in QCD. The t dependence is 
taken from data fits to hard diffractive processes. 

Using the fact that F 2 (x,Q 2 ) ~ ^isiti^Q. one can rewrite Eq. ( |5.34| ) into its 
final form: 

R, - AR2Q 2 B F ^Q 2 )( 1 + v 2 )- ( 5 - 36 ) 

where rf ~ 0.09 — 0.27 for the given Q 2 range. We computed i? 7 , the fractional 
number of events given by Eq. ( |5.36| ), for x between 10~ 4 and 10~ 2 and for a Q 2 of 



2, 3.5, 12 and 45 GeV where the numbers for F 2 were taken from ]75| . Based on our 



analysis of the previous section we use Eq. ( |5.31|) for Q 2 = 2GeV 2 , assuming that 



for Q 2 = 3.5GeV 2 the slope drops by about 1 4- 2 units as compared to Eq. ( [5.31 
to account for the decrease of the transverse size of the gg-pair; for larger Q 2 we use 
Eq. 

We find ~ 1.1 X 10" 3 , 9.9 x 10~ 4 at x = 10" 4 , 10" 3 and Q 2 = 2GeV 2 ; 
i? 7 ~ 1.07 x 10~ 3 , 9.3 x 10~ 4 at x = 10~ 4 , 10" 3 and Q 2 = 3.5GeV 2 ; i? 7 ~ 4.5 x 
10" 4 , 3.78 x 10- 4 2.5 x 10" 4 at x = 10" 4 , 10" 3 , 10" 2 and Q 2 = 12GeV 2 ; and finally 
i? 7 ~ 1.49 x 10" 4 , 1.04 x 10" 4 at x = 10" 3 , 10" 2 and Q 2 = 45GeV 2 . As is to be 
expected, the number of events rises at small x since the differential cross section 
is proportional to the square of the gluon distribution and the total cross section is 
just proportional to the gluon distribution, i.e. , the ratio in Eq. ( |5.34j ) is expected 
to be proportional to the gluon distribution and this assumption is born out by our 
calculation and falls with increasing Q 2 since F 2 does not grow as fast with energy. 
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5.6. The complete cross section of exclusive photon production 

In order to study whether the Bethe-Heitler or DVCS process will be dominant in 
real photon production we need the expressions for the differential cross sections 
first. 

We find that the differential cross section for DVCS can be simply expressed 
through the DIS differential cross section by multiplying the DIS differential cross 
section by i? 7 (see Eq. ( |5.36|) ) which was calculated in the previous section. One can 



see this by observing how F 2 is related to <J to t{l*P) as given in Sec. ^75] and <J to t{l*p) 
to udvcs via Rj in the same section. We then find using Eq. (|5.36|) for i? 7 

<^ = -g?w-yr*- m w*™+tf (5.37) 

with a dvcs = dar> Z as \t=o x 4 using the same exponential t dependence as in the 
previous section and R being the ratio of the imaginary parts of the DIS to DVCS 
amplitudes as computed earlier. 

In writing Eq. ( |5.37| ) we neglected F L (x, Q 2 ) - the experimentally observed con- 
servation of s channel helicities justifies this approximation- so that F 2 — 2xF\. 
y = l—E'/E where E' is the energy of the electron in the final state and r = <f)N+4>e, 
where <pN is the azimuthal angle between the plane defined by 7* and the final state 
proton and the x — z plane and 4> e is the azimuthal angle between the plane defined 



by the initial and final state electron and the x — z plane (see Fig. |5.3|) . Thus <p r is 
nothing but the angle between the 7* — p' and the electrons scattering plane. 

In the case of the Bethe-Heitler process, we find the differential cross section at 
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Figure 5.3: The azimuthal final proton and electron angle in the transverse scattering 
plane. 

small t to be 

da B H __ a 3 sy 2 (l + (l-y) 2 ) 
dxdyd\t\d(p r 7r<5 4 |t|(l — y) 

with t = \t\/Am 2 N and s being the invariant energy. GE{t) and are the electric 

and nucleon form factors and we describe them using the dipole fit 

G E (t) ~ G D (t) = (1 + ^-y 2 and G M (t) = fi P G D (t), (5.39) 

where fi p = 2.7 is the proton magnetic moment. We make the standard assumption 
that the spin flip term is small in the strong amplitude for small t. 

In order to write down the complete total cross section of exclusive photon 
production we need the interference term between DVCS and Bethe-Heitler. Note 
that in the case of the interference term one does not have a spinflip in the Bethe- 
Heitler amplitude, i.e. , one only has F 1 (t), as compared to Eq. (|5.38|) containing a 
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G 2 E (t) + rG 2 M (t) 
1 + r 



(5.38) 



spinflip part, i.e. , Fz(t). The appropriate combination of G\e(£) and Gm(<) which 
yields Fi(t) is 

\G E (t) + ^G M {t) 

(5.40) 



1 + 



1*1 

■liiii 



We then find for the interference term of the differential cross section, where the 
+ sign corresponds to electron scattering of a proton and the - sign corresponds to 
the positron 



aa DVCS+BH 

dxdyd\t\d(p r 



r]a 3 sy(l + (1 - y) 2 )cos(4> r )e- B W 2 F 2 (x, Q 2 



x 



2Qy|t| v /T^R 

G E (t) + TG M (t) ' 
1 + r 



(5.4i; 



The total cross section is then just the sum of Eq. ( |5.37| ),( |5T38| ),( |5.41| ) 



5.6.1. t-dependence of Bethe-Heitler as compared to DVCS for different 

Q 2 ' 

At this point it is important to determine how large the Bethe-Heitler background 
is as compared to DVCS for HERA kinematics, hence, in the following discussion, 
we will estimate the ratio D allowing a background comparison: 



_ < da DV cs+BH/dxdydt > _ 
< d(JBH I dxdydt > 



(5.42) 



with < ... >= Jq W d(j) r . Using the expressions from Sec. |5.6| we compute D and find 



that D > 1 ( See Fig. |5.4| a, |5.6| a) for relatively small y and 0.1 < t < 0.6 with the 
given values of x and Q 2 considered. Note, however, that this does not mean that 
the case for DVCS is hopeless. As it turns out, it is rather advantageous to have 
D < 1 when looking at the interference term which we will do next. 
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It is convenient to illustrate the magnitude of the intereference term in the total 
cross section by considering the asymmetry for proton and either electron or positron 
to be in the same and opposite hemispheres ( we omit the rather cumbersome explicit 
expression but the reader can easily deduce it from Eq. ( |5.37| ),( |5738|) ,( |5.4i|) . ) 



^ _ f-(?/2 d(j) r da D vcs+BH - $l% 2 d(f)rda D vcs+BH 

Jo d(j)rd<JDVCS+BH 

in other words, one is counting the number of events in the upper hemisphere of the 
detector minus the number of events in the lower half, normalized to the total cross 



section. Fig. |5.5| a,b and |5.7| a,b show A for the same kinematics as above and we find 
that the asymmetry is fairly sizeable already for small t and is strongly dependent 
on the energy. Due to this fairly large asymmetry, one has a first chance to access 
nondiagonal parton distributions. We will discuss A in more detail, in particular its 
energy dependence, in the next chapter. 

Note, there is an increased experimental difficulty to measure DVCS if the recoil 
proton is not detected in other words if t is not directly measured. However there 
is a simple, practical way around this problem which we will discuss next. 
5.6.2. DVCS alternative to tagged proton in the final state 

Another interesting process, which can be studied in the context of DVCS, is the one 
where the nucleon dissociates into mass "X" - 7* + p — * 7 + X. Perturbative QCD 
is applicable in this well. In particular the following factorization relation 

should be valid at sufficiently large Q 2 : 



(5.44) 
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The big advantage of the dissociation process as compared to the process where the 
target proton stays intact is that the Bethe-Heitler process is strongly suppressed 
for inelastic diffraction at small t due to the conservation of the electro-magnetic 
current, hence the amplitude is multiplied by an additional factor yj\t\ which is 
basically for the Bethe-Heitler process. Thus, the masking of the strong amplitude 
of photoproduction is small in this case. Since there is already data available on 
J /if} production, this quantity can give us information on how different the slopes 
for the production of massless to massive vector particles are, providing us with 
more understanding on how different or similar the exact production mechanisms 
are. Note that the ratio of the total dissociative to elastic cross section of p meson 
production is found to be about 0.65 at large Q 2 |76|] which is basically of 0(1). 



The same should hold true for J/ip production and in fact this ratio should be 
a universal quantity. This is due to the fact that one has complete factorization, 
hence the hard part plus vector meson is essentially a point and thus for the soft 
part, is does not matter what kind of vector particle is produced. The above said 
implies for Eq. ( |5.44| ) that it also should be of order unity, implying that the order 



of magnitude of the fractional number of events for real photon production to DIS 
remains unchanged even though the actual number of i? 7 might decrease by as much 
as 35%. 
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5.7. Conclusions 

In the above said we have shown that pQCD is applicable to exclusive photopro- 
duction by showing that the ratio of the imaginary parts of the amplitudes of DIS 
to a real photon is calculable in pQCD after specifying initial conditions since the 
derivative in energy of the hard scattering amplitudes can be unambiguously cal- 
culated in pQCD and all the non-perturbative physics can then be absorbed into a 
parton distribution. We wrote down an evolution equation for the imaginary part of 
the amplitude, which can be generalized to the complete amplitude at small x, and 
solved for the imaginary part of the amplitude. We also found that the imaginary 
part of the amplitude for the production of a real photon is larger than the one 
in the case of DIS in a broad range of Q 2 for the reasons as discussed above. We 
also found the same to be true for the full amplitude at small x. We also make ex- 
perimentally testable predictions for the number of real photon events and suggest 
that the number of events are small but not too small such that after improving 
the statistics on existing or soon to be taken data, it would be feasible to access 
the nondiagonal gluon distribution at small x from this clean process. Finally, we 
demonstrated that measuring the asymmetry A at HERA, which is fairly sizable in 
the kinematics in question, would allow one to determine the real part of the DVCS 
amplitude, in other words gain a first experimental insight into nondiagonal parton 
distributions, despite D < 1. 
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Figure 5.4: a) D is plotted versus -t for x = 10" 4 , 10" 3 , 10" 2 , Q 2 = 12 GeV 2 , 
B = 5 GeV~ 2 and y = 0.4. The solid curve is for x = 10~ 4 , the dotted one for 
x = 1CT 2 and the dashed one for x = 1CT 3 . b) D is plotted versus y for the same 
x,Q 2 ,B and -t = 0.1 GeV 2 
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Figure 5.5: a) The asymmetry A is plotted versus —t for x = 1CT 4 (solid curve), 
x = 1CT 2 (dotted curve) and x = 1CT 3 (dashed curve) again for Q 2 = 12 GeV 2 , 
B = 5 GeV~ 2 and y = 0.4. b) A is plotted versus y for the same x,Q 2 ,B and 
-t = 0.1 GeV 2 . 
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Figure 5.6: a) D is plotted versus —t for x = 1CT 4 and 1CT 3 , Q 2 = 3.5 GeV 2 , 
B — 8 GeV~ 2 and y = 0.4. The solid curve is for x = 10" 4 , the dashed one for 
x = 1CT 3 . b) D is plotted versus y for the same x, Q 2 , B and — t — 0.1 GeV 2 
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A plotted versus -t for y = 0.4, 
Q 2 = 3.5 GeV 2 , x=10~* and x=10" : 
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Figure 5.7: a) The asymmetry A is plotted versus —t for x = 10 4 (solid curve), 
x = 1(T 3 (dashed curve) again for Q 2 = 3.5 GeV 2 , B = 8 GeV" 2 and y = 0.4. b) A 
is plotted versus y for the same x, Q 2 , B and —t = 0.1 GeV 2 . 
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CHAPTER 6 
DVCS in DIS at HERA - A Probe of Asymptotia 



6.1. Introduction 



In this chapter based on Ref . |8T| , we will attempt to answer one of the pressing ques- 
tions of high-energy QCD: The energy dependence of the strong amplitude in the 
situation when one or both colliding systems have small size. It is generally agreed 
that the x-range currently available at HERA is not sufficient to test the current 
ideas about the onset of asymptotia via measurements of the parton densities. 

At the same time, the experience in studies of soft processes tells us that the real 
part of the zero angle scattering amplitude, rj, provides us, through the dispersion 
representation over the invariant energy of the collision, with information about the 
energy dependence of the cross section well beyond the energy where n is measured. 
The reason for this is that rj, the ratio of the real to imaginary part of the amplitude 
essentially measures the Ins derivative of the cross section [ |32|j : 



7r d\n(F 2 (x,Q 2 )) 

v =2 dum ■ ( } 

One can also use derivative analyticity relations to derive a more accurate formula 
P2| leading to 



n = tan — . (6.2) 
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for F 2 (x, Q 2 ) oc x~ a . 

In this chapter we propose a new methodology of investigating the energy de- 
pendence of high-energy processes through their real parts and also the shapes of 
nondiagonal parton distributions. We use DVCS as an example which offers us a 
direct way to study nondiagonal parton distributions, though their actual extraction 
from the data is not possible in DVCS due to the fact that the parton distributions 
depend on y\ and y 2 = Vx — x which are dependent variables rather than indepen- 
dent as one would need and thus the inverse Mellin transform of the factorization 
formula cannot be foundQ. The major new result of this chapter is that the current 
successful fits to the F 2 n(x,Q 2 ) HERA data lead to qualitatively different predic- 
tions for the asymmetry, reflecting different underlying assumptions of the fits about 
the behavior of parton densities at x below the HERA range. 



A recent analysis in Ref. [[55|] and the previous chapter have shown that DVCS 
studies at HERA are feasible and we made predictions for the expected DVCS 
counting rate compared to DIS as well as the asymmetry A in the combined DVCS 



and Bethe-Heitler cross section for recent HI data [ITS 



The chapter is structured as follows. In Sec. |6.2| we review the necessary formulas 



of Ref. f5l| and the previous chapter for our analysis. In this context, the formula 



pertaining to the ratio of real to imaginary part of a scattering amplitude at small 
x is of particular importance. We then present the different fits to F 2 (x, Q 2 ) in Sec. 
73] and present the different results for the asymmetry A with respect to t and y, 



1 As was pointed out before this is not true for di-muon production since there, we have two 
independent variables x and £1 . 
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at fixed y and t respectively. Sec. |6.5| contains our conclusions and outlook. 
6.2. Relations between DVCS and DIS 

In order to compute the asymmetry A, we need the ratio of the imaginary part of 
the DIS amplitude to the imaginary part of the DVCS amplitude and the relative 
DVCS counting rate R 7 , expected at HERA in the interesting kinematic regime of 
10" 4 < x < 10" 2 and moderate Q 2 , i.e. , 3.5 GeV 2 < Q 2 < 45 GeV 2 . The relative 
counting rate R 1 is given by 



^^4i^B F2(x ' g2)(1 + r/2) - (0) 

where R is the ratio of the imaginary parts of the DIS to DVCS amplitude as given 
5"5|p| , B is the slope of the t dependence (for more details see Ref. J5S1 and the 



m 



previous chapter.) and r] is the ratio of real to imaginary part of the DIS amplitude, 
i.e. , F 2 (x,Q 2 ), given by Eq. flOl) . 
The asymmetry A is given by |55[ 



. Jli 2 / 2 d(j) r da D vcs+BH - J^j{ 2 d(j) r da D vcs+BH 
A = 1 (.0.4J 

J Q d(j) r d<JDVCS+BH 



where daovcs+BH is given by the sum of Eq. ( |5.37|) , fl5T38| ),( |5.4lD . As explained in 



p5| this azimuthal angle asymmetry is due to the sign change of the interference 
term in the combined DVCS and Bethe-Heitler cross section, when the angle r is 
integrated over the upper hemisphere of the detector as compared to the integration 
over the lower hemisphere of the detector. Although the absolute value of the 



2 We will use the results for R from pa] in our present analysis. 
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parton distributions cannot be extracted from DVCS, the shape of the distributions 
is nevertheless accessible since the real part of the DVCS amplitude is isolated 
through this asymmetry. Therefore, we investigate the influence of different F 2 fits 
on the asymmetry through the relative counting rate which is directly sensitive to 
the ratio of real to imaginary parts of F 2 as shown in Eq. ( |6.3| ) . 

6.3. The different fits to F 2 (x,Q 2 ) 



In the calculation of the asymmetry A we use the recent HI data from Ref. |75| as 



previously used in Ref. [55], a logarithmic fit by Buchmiiller et. al (BH) |f77| , the 



ALLM97 fit [7^] and a leading order BFKL-fit |79j for illustrative purposes 



In the HI data, F 2 behaves for small x as x~ x and hence r] is just |A where 
rj 2 = 0.09 — 0.27 in the Q 2 range given in the previous section. Note that rj has no 
x dependence, for small enough x, and thus depends only on Q 2 . This is not true 
for all of the other fits. 

F 2 in the BH fit takes on the following form 

O 2 074 

F 2 (x,Q 2 ) = 0.078 + 0.364 log(— ^— 2 ) log(- ), (6.5) 

0.5 GeV x 



and hence we find for rj 



,, , Q 2 1/f ,, o2 , 



7] = -0.364 log( o5 ^ eV2 )/F 2 (x,Q^). (6.6) 

Note that this rj has not only the usual Q 2 dependence but depends rather strongly 
on x also, which is not seen in the data for the slope of F 2 . 
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In the ALLM97 fit F 2 at small x takes on the following form 

_Q?_ l 

Q 2 + ml 



Hx, Q 2 ) = 7 J^(F?(x, Q 2 ) + F 2 R (x, Q 2 )), 



with 



where 



F 2 p (x,Q 2 ) = c P (t)x a / {t) 
F 2 R (x,Q 2 ) = c R (t)x R « {t \ 



V 



MS) 



and 



CR (t) = 0.8017 + 0.97307t 3 - 4942 , 

a R (t) = 0.584 + 0.37888t 2 - 6063 , 

c P (t) = 0.28067 + 0.05776 ( 1 + ^ 219re -l), 

a P (t) = -0.0808 + 0.36732 ^* ^ -l). 



xp and x R are given at small x by 



2 



Xp = x(l + 
Xr = X. 



VTLp 

~Q 2 



A 2 = 0.06527 GeV 2 , nip = 49.457 GeV 2 , Q 2 = 0.46017 GeV 2 
0.31985 GeV 2 . rj is then given by 



V = 



ap i aft 
7T apCpXp + CLrCrX r 

"2 c P a£ p + erf 
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4JV e ln(2)a s 

In the case of BFKL where F 2 — x * we find 77 to be 



7T 4iV c ln(2K 
ij = - (6.13) 

Z 71 



with an a s in leading order of 



47T 

* = wzm- (6 - 14) 



6.4. Results for the asymmetry A 



In Fig. |6.1| - |6.3| , we plot the asymmetry A as a function of t and ?/ for fixed Q 2 = 
12 GeV 2 , fixed y = 0.4 and -t = 0.1 GeV 2 and ac = 10" 4 , 10" 3 , 10~ 2 . The slope 
S of the t-dependence for DVCS was taken to be B = 5 GeV -2 whereas for the 
Bethe-Heitler cross section we used the nucleon form factor as used in chapter 5. 
The counting rate _R 7 was appropriately adjusted for the different fits according to 
Eq. ( p.'S\ ). The solid curves in Fig. O - O are our benchmarks^. 



Comparing the BH fit (dotted curves), against our benchmarks we find a strong 
x dependence of the asymmetry in the BH fit as well as different shapes and absolute 
values. The strong x dependence of the BH fit in the ratio of the real to imaginary 
part of F% will make it easy to distinguish this logarithmic fit from a power law fit 
which yields an x independent r]. 

As far as the ALLM97 fit is concerned (short-dash curves), there is hardly a 
difference, as compared to the HI fit in the asymmetry as a function of t and y 



3 Though actual HI data is used, we are still dealing with a leading order approximation and a 
particular model for the nondiagonal parton distributions at the normalization point was used in 
computing i? 7 (see 55 for more details on the type of model ansatz and approximations used.). 
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in absolute value, shape and x dependence, except for x = 10 2 but this is due to 
the approximations we made for xp and Xr which are not that good anymore at 
x = lO- 2 . 

If one compares the BFKL fit (dash-dot curves) to the HI fit one sees immediately 
that the BFKL fit is totally off in almost all aspects and was only included here as 
an illustrative example, not as a serious fit. 

6.5. Conclusions 

In the above we have shown the sensitivity of the exclusive DVCS asymmetry A to 
different F 2 fits and made comments on the viability of each fit. Note that even a fit 
which reproduces F 2 data, as well as its slope, in a satisfactory manner can be shown 
to lead to differences in the asymmetry shape. The sensitivity of the asymmetry 
to y and t will allow us, once experimentally determined, to make a shape fit and 
hence make a shape fit to nondiagonal parton distributions for the first time. 
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Figure 6.1: HI fit (solid curve), the BH fit (dotted curve), ALLM97 fit (short-dash 
curve) and BFKL fit (dash-dot curve) for x = 10~ 4 . a) Asymmetry A versus t for 
fixed y = 0.4. b) Asymmetry A versus y for fixed — t — 0.1 GeV 2 . 
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Figure 6.2: HI fit (solid curve), the BH fit (dotted curve), ALLM97 fit (short-dash 
curve) and BFKL fit (dash-dot curve) for x = 10~ 3 . a) Asymmetry A versus t for 
fixed y = 0.4. b) Asymmetry A versus y for fixed — t — 0.1 GeV 2 
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Figure 6.3: HI fit (solid curve), the BH fit (dotted), ALLM97 fit (short-dash curve) 
and BFKL fit (dash-dot curve) for x = 10 -2 . a) Asymmetry A versus t for fixed 
y = 0.4. b) Asymmetry A versus y for fixed — t — 0.1 GeV 2 
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Appendix A 

Direct Estimation of Sizes of Higher-Order Graphs 



A.l. Introduction 



In this appendix, based on Ref. ||82|| , we want to give concrete examples of distribu- 
tional methods which were employed in an abstract way in chapter 4. The starting 
point for this appendix is formed by the following observations: 

• The only (known) systematic method for calculating scattering in QCD is 
perturbation theory. (Lattice Monte-Carlo methods work in Euclidean space- 
time, and are excellent for calculating static quantities such as masses from 
first principles. But they are essentially useless when a calculation in real 
Minkowski space-time is needed.) 

• In field theory, calculations beyond low orders of perturbation theory are com- 
putationally complex, both because the calculations of individual graphs are 
hard and because there are many different graphs. 

• Hence it is important to make the most efficient use of low-order calculations. 

Since the coupling in practical calculations is not very weak, the accuracy of pre- 
dictions can be ruined by uncalculated higher-order terms. It follows that there is a 
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need to estimate the sizes of the errors. For this one wants quick estimates of terms 
in perturbation theory. The computational complexity of the estimates should in- 
crease as little as possible with the size of the graphs. Indeed, our aim is that one 
only calculates integrals of the form 

U dxx n ln p x. (A.l) 



With suitable methods: 

• One can determine good values for renormalization and factorization scales, 
by asking how to minimize the error estimates. 

• When the estimates get substantially larger than some appropriate "natural" 
size, one would get a diagnosis of a need for resummation of classes of higher- 
order corrections. The diagnosis would include an explanation of the large 
terms and thus indicate the physics associated with the resummation. 

In this appendix, we explain how to start such a program. It builds on work 



first reported in Ref. ||80|| . Our methods treat properties of the integrands of Feyn- 
man graphs, and are therefore directly sensitive to the physics of the process being 
discussed. Some other treatments of these issues discuss the problems in terms of 
the mathematics of series expansions in general, without asking what is causing the 
graphs to be the sizes they are. A particular exception is the work of Brodsky, 
Lepage, and Mackenzie [^3|[|. They use heavy quark loops to probe the actual mo- 
mentum scales that dominate in a particular calculation; this is then used to provide 



1 See also the more recent work of Brodsky and Lu Q and of Neubert |8E 
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a suitable value for the renormalization/factorization scale. But we believe that our 
methods provide a more direct route to answering the question of why the scales 
are what they are and why a calculation gives a particular order of magnitude. The 
issues addressed by methods involving the Borel transform and Pade summation 
address complementary issues |86|j . 

Now, most cross sections in QCD cannot be directly computed by perturbation 
theory; this can only be used to compute the short- distance coefficients that appear 
in the factorization theorem. So we will also treat the specific problems that arise 
in estimating the sizes of the coefficient functions. These functions have the form 
of a sum over Feynman graphs (typically massless), with subtractions to cancel 
some infra-red (IR) divergences. Remaining IR divergences are cancelled between 
different graphs or between different final-state cuts. Two characteristic features 
appear. First, we can obtain estimates for sums of particular sets of graphs, but not 
for the individual graphs, which are divergent. Secondly, the coefficient functions 
are not in fact genuine functions. They are normally singular generalized functions 
(or distributions) and an estimate can only be made for the integral of a coefficient 
function with a smooth test function. 

We show how to estimate the sizes of graphs by a direct examination of the in- 
tegrands. An important part of our technique is an implementation of subtractions 
directly in the integrands, both for the subtractions that implement counterterms 
for ultra-violet (UV) renormalization and for the IR subtractions that are used in 
short-distance coefficient functions. In order to explain our ideas, we will examine 
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P-k 



Figure A.l: One-loop self-energy graph. 

two examples: (1) the one-loop self-energy graph in (0 3 )4 theory, and (2) a particular 
set of graphs for the Wilson coefficient for deep-inelastic scattering. Our estimates 
are in the form of approximations to ordinary integrals that are absolutely conver- 
gent. This is in contrast to the original integrals, which are typically divergent in 
the absence of a regulator. Thus a by-product of our work will be algorithms for 
computing graphs numerically in Minkowski space-time, which may have relevance 



to work such as Ref. [87j. As an illustration of how estimations can be carried out, 
even analytically, for a measurable quantity, we will estimate the size of the Wilson 
coefficient for the structure functions Ft and Fl in the last section. 

A. 2. Euclidean self-energy in (ft 3 theory 

In this section, we give a representation of the one-loop self-energy graph of Fig. |A.l 



in 3 theory in four dimensions. A particular renormalization scheme is used, which 
we relate to ordinary MS renormalization. Then we show how to estimate the size 
of the graph from elementary integrals and hence how to choose the renormalization 
scale suitably. 
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A. 2.1. Renormalization 

The only complication in computing the graph of Fig. |A.1| is its UV divergence, 



given that we choose to work in Euclidean space. We give a representation of the 
graph that is an absolutely convergent integral over the loop momentum itself in 
four dimensions and in which the renormalization is explicitly "minimal" . This last 
term means that the counterterm for the (logarithmic) divergence is independent of 
the mass and the external momentum. It is a very useful property when one wants 
to take zero-mass limits, etc. For a general divergence, the counterterm would be 
polynomial in masses and external momenta. 
Our representation of the graph is: 

KP > 2(2vr) 4 7 \{k 2 + m 2 ) [(p - k) 2 + m 2 } k 4 J ' 1 

We recognize in the integrand a term that is given by the usual Feynman rules 
(in Euclidean space-time), and a subtraction term. The subtraction term is the 
negative of the asymptote of the first term as k — > oo, so we term our procedure?] 
"renormalization by subtraction of the asymptote" . A cut-off is applied to prevent 
the subtraction term giving an IR divergence at k = 0; the cut-off does not affect 
the k — > oo behavior and therefore does not affect the fact that the UV divergence 
is cancelled. To see that Eq. ( |A2|) is equivalent to standard renormalization, one 



simply applies a UV regulator, after which each term can be integrated separately. 

The first term is the unrenormalized graph and the second term is a p-independent 

2 See |8Cj for a previous account. A formalization of such ideas (to all orders of perturbation 
theory) was given earlier by Ilyin, Imashev and Slavnov |8S| , and later by Kuznetsov and Tkachov 
& 
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counterterm. 

Evidently, the integral is absolutely convergent, and can therefore be computed 
by any appropriate numerical method. (It can also be evaluated analytically. But 
this is not interesting to us, since we wish to obtain methods that work for integrals 
that are too complicated for purely analytic methods to be convenient or useful.) 

The counterterm is in fact the most general one that is independent of m and p, 
since any other renormalization counterterm can differ only by a finite term added 
to the integral that is independent of m and p, and a change of the cut-off p c is 
equivalent to adding such a term. We can relate the counterterm to the commonly 
used MS one simply by computing the counterterm alone, with dimensional regu- 
larization: 



The result is that setting p c equal to the scale p of the MS scheme gives exactly 
MS renormalization. In general, we would find that p c would be a factor times p, 
or equivalently that we should set p c = p and then add a specific finite counterterm 
to the graph. 

Notice that the integral to relate our renormalization scheme to the MS scheme 
is algorithmically simpler to compute analytically than the original integral. There 
is always the possibility of adding finite counterterms. Moreover, the precise form of 
the cut-off is irrelevant to the general principles. One can, for example, change the 
sharp cut-off function 6{k > p c ) to a smooth function f(k/p c ) that obeys /(oo) = 1 
and /(0) = 0. Such a function would probably be better in numerical integration. 




(A.3) 
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Of course, our method as stated is specific to one-loop graphs. But it is an idea 



that has been generalized |88|, |89| to higher orders. 
A. 2. 2. Estimate 



We now show how to estimate the size of the integral Eq. (|A.2 ). To give ourselves 
a definite case, let us choose p 2 < m 2 . We obtain the estimate as the sum of 
contributions from k < m and from k > m. Since the renormalization counterterm 
is designed to subtract the k —>■ oo behavior of the unrenormalized integrand, we 
regard it as a 5-function at infinity and therefore to be associated completely with 
the k > m term in our estimate. 

In the region k < m, our estimate is obtained by replacing each propagator by 
1/m 2 so that 

Contribution from k < m ~ / d 4 k — - 

32vr 4 Jk<m m 4 

- & < A ' 4 > 

This factor is the product of g 2 /327r 4 for the prefactor and tt 2 /2 for the volume 
of a unit 4-sphere. As advertised, we have had to calculate no integral that is 
more complicated than a simple power of k. The approximation of replacing the 
propagators by 1/m 2 leads us to an over-estimate of the integral, but not by a great 
factor, since we are in a region of small momentum. 

The estimate for k > m is obtained by replacing the propagators by their large 
k asymptote: 

g 2 [ j4 ,[l 6{k>fi c y 
k 4 k 4 



Contribution from k > m ~ I d k 

327T 4 Jk>m 
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-Ztt 



32vr 4 



00 dk r°° dk 

1 k J a h 



(A.5) 



Since we are taking the difference of two terms, we must be careful about the errors, 
which are of order 



To understand the structure of the result, let us examine how the original integral 
( |A.2| ) appears after integrating over the angle of k: 

Ifb/o ^[ A (k,P,m)-9(k>i, c )]. (A.7) 

The function A is the angular average of k 4 times the two propagators. It approaches 
as k — > 0, so that the integral is convergent there, and it approaches unity as 
k — > 00, which would give a UV divergence were it not for the subtraction. 

We can represent this situation by the graphs of Fig. [A.2|. The transition region 



for the unsubtracted integrand, where it changes from being to 1, is around k = m, 
to within a factor of 2 or so; this is evident by examining the denominators. By 
setting fi c to be m to within a factor of 2, we achieve the following: 

• The integrand is less than unity everywhere. 

• It is concentrated in a shell of thickness of order m around k = m. 
Thus we can say that the natural size of the graph is the product of 

• The prefactor c/ 2 /32tt 4 . 

• The surface area of a unit 4-sphere: 2tt 2 . 
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(a) 



Ink 




Figure A. 2: (a) Scaled integrand for unrenormalized self-energy graph, (b) Subtrac- 
tion term, (c) Total if fj, ~ m. (d) Total if ji ^> m. 
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The width of the important region, about unity, in units of m. 
A factor of m to the dimension of the integral, i.e., 1. 

A factor less than unity, say 1/2, to allow for the fact that, after subtraction, 
the integrand in Eq. (A. 7) is smaller than unity and that there is a cancelation 



between negative and positive pieces. 



That is, the natural size is g 2 /32ii 2 , if /i is reasonably close to m. 



If \x is not close to m, then we get a long plateau in the integrand — see Fig. |A.2 . 
The height of the plateau is unity, and its length is ln(/A/m) to about ±1 in units 
of In k. This clearly gives a larger-than-necessary size, and an optimal choice of the 
MS scale is around m. 

One should not expect to get an exact value for the scale fi. A physical quantity 
in the exact theory is independent of fi, and any finite-order calculation differs from 
the correct value by an amount whose precise value is necessarily unknown until 
one has done a more accurate calculation. If one is able to estimate the size of the 
error, as we are proposing, then an appropriate value of fi is one that minimizes the 
error. Given the intrinsic imprecision of an error estimate, there is a corresponding 
imprecision in the determination of \x. One can expect to identify, without much 
work, an appropriate scale /i to within a factor 2, and, with a bit more work, to 
within perhaps 50%. These estimates just come from asking where the transition 
region in the integrand is, and by then obtaining an answer by simple examination 
of the integrand. But one cannot enter into a religious argument of the wrong kind 
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as to whether the correct scale is 1.23m as opposed to 1.24m, for example.^ By 
definition an error estimate is approximate. 

Since we have not yet investigated how to estimate even higher-order graphs, we 
are making the reasonable conjecture that the properties of graphs do not change 
rapidly with order. Then our estimate that /i should be close to m will ensure that 
higher-order graphs are of the order of their natural size. 

A. 2. 3. Implication for QCD 

The same arguments applied to similar graphs in QCD show that the natural ex- 
pansion parameter in QCD is 

OL 

— - x (group theory) x (factor for multiplicity of graphs). (A. 8) 

4:71 

These arguments rely on being able to show that in the dominant part of the range of 
integration all lines have approximately a particular virtuality and that the relevant 
range of integration is a corresponding volume of momentum space. 

In the general case we cannot expect to get a much smaller result, but we can 
expect that in some situations the properties of the dominant integration region(s) 
will not be so good. So what we need to do next is to analyze more interesting 
graphs in QCD. This we will do in the next section. 

In general, when we get corrections in QCD that are substantially larger than the 

natural size given above, it must be either because the integrand is excessively large, 

3 Brodsky and Lu obtain very precise estimates of a suitable scale. Their rationale is the 
elimination of IR renormalons in the relations between IR-safe observables. This is a concern with 
very high-order perturbation theory, an issue that we do not address. 
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or because there is no single natural scale, or because we have not chosen a good 
scale.f] If the integrand is especially large relative to the natural unit, or if there is 
no single natural scale, then we should investigate in more detail the reasons, and 
derive something like a resummation of higher-order corrections |K|. In precisely 
such situations, one does indeed have to compute high-order graphs. At the same 
time, there is no need to compute the complete graphs in all their gory detail, but 
only their simple parts. 

A. 3. Wilson coefficient for deep-inelastic scattering 

The factorization formula for the leading-twist part of a deep-inelastic structure 
function F{x) is 

Fix) = J jf(t)H*/t)- (A-9) 

Here, /(£) is a parton density, and F is the short-distance coefficient ("Wilson 
coefficient"). We have suppressed the indices for the different structure functions 
(Fx, F 2 , etc.) and for the parton flavor. The coefficient function F(x/£) is obtained 
from Feynman graphs for scattering on a parton target with momentum £p, where 
p is the momentum of the hadron target. Subtractions for initial-state collinear 

singularities are applied to the Wilson coefficient and the massless limit is taken. 

4 Our use of the word "natural" may suggest that we are proposing to estimate higher-order 
corrections simply by multiplying the appropriate power of the natural expansion parameter by 
the number of graphs. This is not what we mean. We are arguing first that the sizes of graphs 
can actually be estimated fairly simply, and that the natural expansion parameter is a useful unit 
for these estimates. Secondly, we show that in the most favorable cases, graphs are less than or 
about unity in these natural units. Finally, we argue, in the next section, that general kinematic 
arguments about the physics of a graph are useful in diagnosing cases where graphs are large in 
natural units. 
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(a) © 



Figure A. 3: The cuts of a one-loop graph for the Wilson coefficient for deep-inelastic 
scattering. 

As an example, we will examine the contribution to the Wilson coefficient from 



the diagrams in Fig. |A.3| . These diagrams are the two possible cuts of a particular 
uncut one-loop graph, and since we will need to use a cancelation of final-state 
interactions, we must consider the sum of the two cut graphs as a single unit. Note: 

• We will apply a subtraction to cancel the effect of initial-state collinear in- 
teractions where the incoming quark splits into a quark-gluon pair which are 
moving almost parallel to the incoming particle .Q 

• There will be soft-gluon interactions and collinear final-state interactions. 
These will cancel after the sum over cuts. (In a more general situation, a 



5 According to the factorization theorem, the subtractions cancel all the sensitivity to small 
momenta, i.e., to the initial-state collinear interactions. The subtractions are of the form of terms 
in the perturbative expansion of the distribution of a parton in a parton convoluted with lower- 
order terms in the coefficient function. (See, for example, ]9l], |34| for details.) Of course, both the 
partonic cross section and the subtraction term have to be properly renormalized. We will call the 
subtraction terms eikonal because of the particular rules involved in their calculation, for graphs 



such as Fig. A. 3 
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sum over a gauge-invariant set of graphs is necessary to get rid of all soft- 
gluon interactions.) 

The Wilson coefficient is a distribution (or generalized function) rather than 
an ordinary function of x/£. Thus it is useful to discuss only the size of the 
coefficient after it is integrated with a test function, but not the size of the 
unintegrated coefficient function. The parton density /(£) provides a ready- 
made test function that has a physical interpretation. 



The first graph of Fig. |A.3|, which has a virtual gluon, has a 4- dimensional 



integral, but the second graph, with a real gluon, has only a 3-dimensional 
integral. Thus the cancelations associated with the sum over cuts can only be 
seen after doing at least a 1-dimensional integral. 

It is useful to visualize the process in space-time, Fig. [A.4|, and to use light-front 



coordinates (+, - T) (defined by V ± = (V° ± V 3 )/V2). O ur axes are such that the 
incoming momenta for the hard scattering are: 

£pf = (& + , 0, Or), <f = (-xp + , T J . (A.10) 



Also, we find it convenient to parameterize the gluon momentum in Fig. |A.3| in terms 



of two longitudinal momentum fractions, u and z, and a transverse momentum k^, 
as follows: 

W = (u£(l-z)p + , zq~, k T ). (A.ll) 

Thus z is exactly the fraction of the total incoming minus component of momentum 
that is carried off by the gluon, while u is a scaled fraction of the plus compo- 
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Figure A. 4: Space-time structure of deep-inelastic scattering, in the center-of-mass 
frame of the virtual photon and the struck quark. The solid line is the almost (light- 
like) world line of the incoming quark. The dashed line is the world line of the single 
struck quark in the lowest-order (Born) graph for the hard scattering. 



139 



nent. The scaling is somewhat unobvious, but it has the effect that positive energy 
constraints on the final state restrict each of u and z to range from to 1. 

Now we summarize how the calculation of the real and virtual graphs and of the 
subtraction graphs ||91||, contributing to the Wilson coefficients, is to be carried out: 



Using the Feynman rules for cut diagrams, we write down the momentum 
integral with the appropriate 5-functions. Then we contract the trace over 
Dirac matrices with the appropriate transverse, longitudinal or asymmetric 
tensor on the photon indices to obtain the structure function Fi, F2, etc. For 
the sake of a definite simple example we choose to contract with —g^v, which 
in fact gives the combination 3Fi — F 2 /2x. 

We express the results of the calculations in terms of the light-cone components 
of k (the internal gluon momentum), of p (the incoming quark momentum), 
and of q (the incoming photon momentum). 

In the graphs with real gluon emission, we use the two (^-functions to perform 
the integrals over ky and over the fractional momentum £ entering from the 
parton density. This leaves a 2-dimensional integral. 

In the virtual graphs, we use the one 5-function to perform the integral over 
the fractional momentum £ entering from the parton density. We then perform 
the kr integral analytically. Again we have a 2-dimensional integral. 

We change the integration variables to the scaled dimensionless variables u 
and z defined in Eq. ( |A.11|) . This gives us an overall factor, just like the 
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Final state 
Singularity 



> u 



Initial State Singularity 



Figure A. 5: Integration region and position of singularities for Fig. A.3| (a). 



g 2 /167r 2 in the self energy, times an integral over roughly the unit square in u 
and z. 



A.3.1. Real gluon 



For the real gluon graph Fig. |A.3| (a) it is well known that the integrand has the 
following singularities in the massless limit: 

• Initial-state collinear singularity on z = 0. 

• Final-state collinear singularity on u = 0. 

• Soft singularity k = at the intersection of the previous two singularities, i.e., 
at u — z — 0. 



(See Fig. |A.5| .) In accordance with the standard recipe for constructing the Wilson 
coefficient, we subtract the initial-state collinear singularity. The subtraction term 
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itself has a UV divergence, which we choose to cancel by using the same method as 
we used for the UV divergence of the self-energy graph. This givesQ 



9 f 

T^ Cf / dz , 
%it z Jo Jo 



l-x 



du f 



X 



l-u 



1 - z 



zu 



[Z < Z r , 



ZU 



(A.12) 



The cut-off z cnt on the subtraction term is analogous to the cut-off fi c we used for 
the UV counterterm for the self-energy graph. The value of z cut needed to reproduce 
the MS prescription can be found by a simple calculation from the Feynman rules 
for parton densities. But we will not use this result here. Rather, we will aim at 
calculating an appropriate value for z cnt to keep the one-loop correction down to a 
"normal size" (and, most importantly, whether it is possible to find such an appro- 
priate value at all). This effectively amounts to a choice of factorization scheme. 
Once a suitable value for z cnt has been obtained, it is a mechanical matter to trans- 
late it to a value for /ly^ (or to a value of the scale /i in any other chosen scheme). 
The calculation may also result in a need for an extra finite counterterm. 

The subtraction in Eq. ( |A.12[ ) has evidently accomplished its purpose of can- 
celling the initial-state singularity. But we are still left with the singularity on the 
line u = 0. This singularity will cancel against a singularity in the virtual graph, as 
we will now see. 

A. 3. 2. Sum of virtual and real graphs 



Next we compute the virtual graph of Fig. |A.3| (b), following the same line as in 

the previous subsection. We will construct an integral in the same variables as the 

6 We have chosen the overall normalization of the graphs to be such that the lowest-order Born 
graph gives just f(x). 



142 



real graph. The reason why we do this is that the cancelation of the divergence 
at u = will be point-by-point in the integrand. This can be seen from the proof 
by Libby and Sterman |Kj . They treat a general case of final-state interactions, of 
which our example is a particular case. They first treat one integration analytically, 
with the aid of the mass-shell conditions for the final state, in such a way that the 
integrations for graphs related by different positions of the final-state cut then have 
the same dimensions. After that, the cancelation between the different graphs is 
point-by-point in the integrand. 

This implies that we need to perform the £ and k T integrals. We do the convo- 
lution with £ by the mass-shell 5-function, which now gives £ = x. Then we do the 
k^-integral analytically. (This is not the most trivial integral, but it works conve- 
niently with our choice of variables.) An example of the type of integral encountered 
is: 

Jo [Q 2 (z - l)u -k 2 T + ie] [Q 2 (u - l)z - k 2 + ie] [Q 2 uz - k 2 T + ie] ' ^ A ' 13 ^ 

In the virtual case, no longer does a positive energy condition restrict the range of 
u and z. Nevertheless, it is convenient to split up the result into a piece inside the 
square < u, z < 1 and a piece from outside the square. It is sufficient to examine 
the contribution within the square. After subtracting the collinear singularity and 
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adding the result from the real graph we getQ 

1 - z 6{z < z cat ) 



Q 2 f 1 f 1 
— -C F / dudz 
8n 2 Jo Jo 



zu zu 

+ contribution from outside < u, z < 1. (A. 14) 

The UV divergence (k — > oo) is from outside the square and is cancelled by sub- 
tractions just as for the self-energy graph. 

We see that the singularity at u = for fixed z has cancelled, even at the 
intersection of the two singular lines. However, how good the cancelation is depends 
on the test function. If f(x) is slowly varying, the cancelation is good over a wide 
range of u. But if f(x) is steeply falling as x increases, the cancelation is good only 
over a narrow range of u, and we are left with an integrand that behaves like 1/u 
for larger u. For fixed z we then have an integrand like that in Fig. [A. 6 . 



This kind of behavior is quite typical when singularities are cancelled between 
different graphs with different final states. The integrand is a non-trivial distribu- 
tion, and how good the cancelation of the singularity is depends on properties of 
the test function. The cancelation only occurs after integration over a range of final 
states. This is in distinct contrast with the case of the initial-state singularity for 
which we have constructed an explicit subtraction. The cancelation of the z — * 
singularity occurs before the integration with the test function. 

Immediately we also get complications in computing the typical size of the graph. 



7 



I n the upper limit on u is 1 — x, but we replace the limit by 1 when we copy the formula 



into ( A. 14 ). This change is innocuous since the limit 1 — x arises from the fact that the parton 



density /(£) is when £ > 1, and hence that f(x/{l — u)) is when u > 1 — x. The limits on u 



that result from positivity of energy of the two final-state lines in Fig. A. 3(a) are just < u < 1. 
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A 



1 — 



Figure A. 6: The integrand of Eq. (|A.14j) as a function of u at fixed z when f(x) is 
a steeply falling function of x. 

In particular, when f(x) is steeply falling, we should therefore expect a size for the 
integral that is much larger than the natural size we defined earlier. 

A. 3. 3. Estimate of size 

First, as a benchmark case, let us assume that f(x/(l — u)) is slowly and smoothly 
varying when u increases from to 1/2, and thatQ z cut is around 1/2. Then the size 
of ( 1A.14j ) can be estimated as 



9 2 , 9 2 
— - x (group theory) x (area of unit square in (u,z)) x fix) ~ — of( x )- (A. 15) 
on 2 on 2 



8 The reason for using 'u — 1/2' in these criteria rather than, say, u = 1 is the same as for 
using fi = m in the calculation of the self-energy. It is a rough attempt to optimize the errors 
without using the details of the integrand, since u = 1/2 is midway between the singularity at 
U = and the approximate edge of the region of integration. The integration region for the real 
graph extends to u = 1 — x, while that for the virtual graph extends beyond u — 1. Analogous 
reasoning applies to z and the relation between the real graph and its eikonal approximation even 
though, in this case, one has to look at the whole integrand since both test functions are of the 
form / {-^ 
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A way of obtaining this result with the same method as we used for the self-energy 
is to write 



Ci? x range of fc, i.e., 27T Q x size of integrand, . (A. 16) 

(2vr) 4 v y V / 

Given that the lowest-order graph is f(x), and that the integrand varies in sign, all 
this implies that the contribution of this graph (with all the cuts and subtractions) 
is probably somewhat smaller than g 2 /6ir 2 times the lowest-order graph. In other 
words, the simplest estimate for the graphs, Eq. ( |A.8|) , is a valid estimate in this 
situation. 

There are a modest number of graphs, so this result would have very nice impli- 
cations for the good behavior of perturbation theory: the real expansion parameter 
of QCD would be a s /n, which is a few per cent in many practical situations. 

Unfortunately, the conclusion is vitiated when f(x) is steep, as is often the case. 
Consider the parton density factor times the 1/u factor, relative to the lowest-order 
factor f(x). In the limit u — > 0, this is 

/ fe) ~ fW xf{x) 
uf{x) f{x) 

This factor should be of order unity, if the previous estimate of the size, Eq. ( |A.15[) , 
is to be valid. 

But if the logarithmic derivative xf'(x)/f(x) is much bigger than 1, then we 
have to change our estimates. Consider a typical ansatz for a parton density: 

f{x) oc (1 - x) 6 , (A.18) 
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for which 

x_f 
f 



(A.19) 

1 — x 



This is 6 when x = 1/2, and goes to infinity as x — > 1. Clearly our estimate in 
Eq. ( |A.15| ) is bad. Moreover the u — > estimate, Eq. ( |A.17p , is only approximately 



valid when f(x) does not change by more than a factor 2 (roughly). Once u gets 
larger than the inverse of the logarithmic derivative, the f(x/(l — u)) term in Eq. 
(|A. 14|) is no longer important, and we get the result pictured in Fig. |A.6| : we have 



basically al/ti form with a cut-off at small u. This is a recipe for a large logarithm, 
with the argument of the logarithm being the large logarithmic derivative. 
To make the estimate, it is convenient to define 

5u = . (A.20) 

2f'{x)x v ; 

From Eq. (|A.17|) , we see that Su is approximately the change in u to make f(x/(l — 
u)) fall by a factor of 2. We interpret 5u as the value of u at which the final-state 
cancelations become "bad" . 

Next we note that for normal parton densities f(x) is a decreasing function of x. 
So a simple useful estimate can be made by making the following approximation: 



x \ „ J f( x 



1 " 



25u 

if u > 2Su 



if // < 2<)u. 



Therefore, in Eq. ( ^Tg ), we can replace $ %f(x/(l - u)) by f* Su ^f(x). 

Our estimate for Eq. ( |A.14j ) is the sum of contributions from the following re- 
gions: 
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• < u < 2Su: The value of the integrand is about —f(x)/26u times a function 
of z. For good choices of z cnt , the function of z is less than about unity, but 
with an indefinite sign. Thus we obtain a contribution of about f(x) in size. 

• 25u < u < 1: The integrand is now approximately —f(x)/u, again times a 
mild function of z. We therefore obtain a contribution of about f(x) \n(l/2Su). 

• Exterior of unit square: Here, the only contribution is from the virtual graph, 
and we have no final-state singularity. Hence the naive estimate of unity is 
valid. (The precise value, when /i — Q, is in fact somewhat larger.) 

All of these are to be multiplied by the prefactor g 2 CV/87r 2 . So we obtain a total 
contribution of 



where each term represents an estimate, valid up to a factor of 2 or so. The con- 
tributions from within the unit square may have either sign, depending on the cut 
in the collinear subtraction, while the contributions from outside the unit square, 
from virtual graphs only, have a negative sign. As an explicit indication that our 
estimates are valid for the sizes but not the signs of the graphs, we have inserted 
a ± sign in front of each term. This estimate assumes that renormalization of the 
UV divergence of the virtual graph is done at the natural scale fi ~ Q, and that 
renormalization of the parton densities is done so that it corresponds to z cut of about 
1/2. It also assumes that f(x) falls steeply enough for 5u to be less than about 1, 
as is typically true. 




(A.22) 
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A. 3. 4. Interpretation 

It is obvious that there is a logarithmic enhancement in Eq. ( |A.22| ) whenever 5u is 
small. Our calculation is, of course, no more than a rederivation of the standard 
observation that there are large logarithms in the x — > 1 region. What our derivation 
adds is to show that it is not so much the limit x — > 1 that is causing the problem as 
the steepness of the parton densities. Moreover we have given a numerical criterion 
for when the correction begins to be larger than what we called the natural size 
for higher-order corrections. In other words we have shown how to estimate the 
constant term that accompanies the logarithm. Moreover this is all presented in the 
context of a general method for obtaining estimates of the sizes of graphs. 

It is perhaps clear that, with sufficient foresight, one could have predicted the 
large corrections merely from the observation that the derivation of the factorization 
theorem requires the cancelation of final-state divergences between different final 
states. 

There is in fact another source of large corrections that will make its effect felt 
in even higher order. This is a mismatch in the scales needed for renormalizing the 
parton densities. We have renormalized these by using a value of z cnt that must be 
about 1/2 to avoid making the contribution of the graph unnecessarily large. In the 
case of the real graph, we can translate this to a scale of transverse momentum by 
using the mass-shell condition for the gluon: 

4 = Q 2 uz(l - z)^ = Q 2 UZ ^~ Z \ (A.23) 
x 1 — u 
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Evidently, whenever small values of u are important, small values of kr (relative to 
Q) will be important. This does not affect our one- loop calculation. But in higher- 
order correction, the virtuality of some internal lines will be controlled by the value 
of kx-i and hence there will be mismatches between the scales needed at different 
steps in the calculation. 

Once one has diagnosed the problem, we see that a proper solution lies in more 
accurately calculating the form of the Wilson coefficient near its singularity. This 
subject goes under the heading of resummation of large corrections |9(J . 

A. 4. More detailed estimation of Ft,l to one-loop order 

In the following we will estimate the sizes of the one-loop corrections to the struc- 
ture functions Ft and Fl without doing actual calculations of Feynman diagrams 
by giving a recipe of how to construct the estimates from general principles and 
kinematic considerations. However, we will present the recipe in the context of an 
actual set of Feynman graphs. By using the calculations of the graphs in Sec. |A.5| , 
we will verify that these "simple-minded" estimates are actually valid. 

One obtains Ft,l from the well-known hadronic tensor W^ u by projecting out 
the "transverse" and "longitudinal" pieces via: 

F T = -9^W^ = 3F x -£ 



Q 2 PnPv^ rlxv _ p , F 2 

which give, for example, 



Fl = ^j^ w , v = _ Fl+ < (A.24) 
p ■ q z 2x 



2 



~9l±V n 

p ■ q A 
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2, 

' " ' W^. (A.25) 



A. 4.1. Estimation of F L 

The recipe for estimating Fp is the following: 

• The singularities we encounter (UV, collinear and soft) are all in the form of 
a factor times the Born graph, and the Born graph has no longitudinal part. 
Therefore the one-loop graphs for Fl have no UV, soft or collinear singularities. 

• Since the parton densities are falling with increasing x, the size of a graph is: 

g 2 

x Cp x f[x) x range of u x range of z. (A. 26) 



• The range of z is 1. 

• The range of u is 5u. 

It is elementary to show that the self-energy and vertex graphs (whether real or 
virtual) give a zero contribution to F L . Therefore, our result for F L is: 

F L = f 1 f{x)5u. (A.27) 

Let us now check whether our intuition has guided us in the right way. We use 
Eq. ( |A.37|) for the contribution to the coefficient function for Fl. We approximate 



the u integral by 

/ duf ~ 6uf ( x ^ ( A ' 28 ) 

which is appropriate for a typical parton density, which falls with increasing x. The 
z integral gives a factor 1/2, and we recover Eq. ( A.27[ ), which we obtained by more 
general arguments. 
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Notice that Fl is generally rather smaller than what we have termed the natural 
size, because of the 8u factor: i.e., because of the restricted phase space available. 
There are no enhancements due to final-state singularities. 

Our estimation methods can be applied in two ways. One is to estimate the 
graphs without having explicit expressions for the graphs; one just searches for the 
possible singularities that would prevent the natural size of a graph from being its 
actual size. The second way of using the methods is to examine the expressions 
for graphs with the knowledge of the singularity and subtraction structure and to 
perform a more direct estimate of the sizes. This is useful since the integrals must 
often be performed numerically, e.g., whenever one convolutes with parton densities 
that are only known numerically. Additional information that is now obtained 
concerns the typical virtualities, etc., of the internal lines of the graphs (compare 
Neubert's work fl85|). This enables a diagnosis to be made of the extent to which a 
problem is a multi-scale problem and therefore in need of resummation. 

A point that we have not addressed is the estimation of the size of the trace 
of a string of gamma matrices, for example. In the standard formula for such a 
trace involves a large number of terms. Nevertheless, it is evident from the above 
calculations that there are cancelations. The final result is that we obtain, relative 
to corresponding numbers for a scalar field theory, a small factor (1 to 4) times a 
standard Lorentz-invariant quantity for the process. 

Evidently, more work on this subject is needed. But the issue of the size of 
the numerator factors from traces, etc., affects completely finite quantities, such as 
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the one-loop coefficient function for Fl, just as much as quantities with divergences 
that are cancelled. However, it is the latter quantities that have the potential for 
especially large corrections. Numerator factors result in magnitudes common to all 
graphs. 

A. 4. 2. Estimation of F T 



We have already calculated one graph, Fig. |A.3j, for Ft, so we only need to sum 



marize the method and apply it to the remaining graphs. The general procedure 



is: 



Since we have singularities in our graphs, they have to be cancelled. In Wilson 
coefficients such as we are calculating, there are explicit subtractions for the 
collinear initial-state singularities. Then there are explicit subtractions for 
UV divergences (both those associated with the interactions and those needed 
to define the parton densities and that therefore enter into the initial-state 
subtractions). Finally there are final-state singularities that cancel between 
real and virtual graphs. 

We consider separately the integrations inside and outside the square < 
u, z < 1. 

A term of the order of the natural size arises from outside the square < u, z < 
1. This comes from purely virtual graphs, with their collinear subtraction. 

Inside the square, a real graph without a final-state singularity contributes an 
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amount of the order of the natural size times Su. This reflects the restriction 



on the range of integration imposed by the parton density. 

Similarly a virtual graph without a final-state singularity contributes a term 
of the natural size. 

Finally, the sum of real and virtual graphs with a final-state singularity 
contributes a term of the order of the natural size enhanced by a factor 
2 + ln(l/2<5u), just as in our estimate of Fig. |A.3| . 



We have the graphs of Fig. A. 3, their Hermitian conjugates, the cut and uncut 



self-energy graph of Fig. |A.7| , below, and the ladder graph, Fig. |A.8| . Since 5u < 1 
typically, the ladder graph gives a small contribution, and it is sufficient to multiply 
the estimate of the Fig. [A.3| by 3: 



2tt : 



:/(*) 



2 + ln 



1 



2Su 



(A.29) 



To see how this compares with the results from the actual graphs, we use those 
in Sec. IA.5. 



The cut and uncut self-energy graphs, Fig. |A.7| below, have final-state singu- 
larities, as can be seen in Eq. (|A.35[) . We found it convenient to use z and kx as 
integration variables. There is a singularity at kq> = in each individual graph. In 
analogy to the definition of 5u, Eq. (|A.20|) , we define 



5k' 



-Q 2 f(x) 
2f'(x) 



(A.30) 
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and by following the same steps as we applied to the ladder graph, we find the 
following estimate for the contribution of Fig. [A.7|: 



q 2 

' C F f(x) 



16tt 



1+J ^ 



25k 



2 



(A.31) 



We have inserted a — sign in this estimate; it is fairly easy to see that the coefficient 
of the logarithm is negative. 

For the ladder graph, whose contribution to Ft is in Eq. ( |A.36| ), there is only 
an initial-state singularity and that is cancelled by an explicit subtraction. If z cut 
is around 1/2, then we can apply the same reasoning as for the longitudinal part of 
the ladder graph, and we find an estimate 

± £jf(x)5u, (A.32) 



where, as in Eq. (|A.22 ), we use the ± to indicate that our estimates do not determine 



the sign of the contribution. 

We have already examined the cut and uncut vertex graphs, in Eq. ( |A.22j ). When 
we multiply this by 2 (to allow for the Hermitian conjugate graphs) and add the 
self-energy and ladder contributions, Eqs. ( |A.31| ) and ( |A.32| ), we get somewhat less 



than our original estimate, Eq. ( |A.29|) , provided the renormalization mass [i is in a 
reasonable range. This lower value is because the self-energy graph is simpler than 
the vertex graph. 

A. 5. The remaining longitudinal and transverse contributions to the Wil- 
son coefficient 
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+ 



(a) 




(b) 



Figure A. 7: (a) Real self-energy graph, (b) virtual self-energy graph. There is a 
second virtual graph that is the conjugate of graph (b). 

A. 5.1. Transverse part of the self-energy 

The real part of the self-energy is computed along the same lines as mentioned in 
the main body of the text, except that we chose to integrate over u instead of kx as 



a matter of convenience. The graph and its possible cuts are shown in Fig. |A.7| . As 
far as the virtual part is concerned, one notes first that the general structure of the 



quark self-energy for zero mass is [92 



where B(0) is computed via: 

B(0) 



B(0)p, 



1 Tr 

4 y d P + ■ 



(A.33) 



(A.34) 



The result for the virtual graph is then given by simply multiplying the Born result 
by -B(O) and convoluting with the test function. The complete result according to 



our prescription is: 



8tt 2 



dz 

JO 



^i^H^f^H- (A - 35) 
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Figure A. 8: The ladder graph for the Wilson coefficient. 

where ^ = 1 + after having made the change of variable k\ — > k\z(\ — z), which 
also gives k\ max = 9Jlzll _ Renormalization by subtraction of the asymptote is used 
in our formula with a cut-off on the momentum integral for the virtual graph of fi 2 . 
As one can see there are no initial-state singularities simply because propagator 
corrections do not induce initial-state singularities as vertex corrections do. The 
final-state singularities are still there; thus one needs to look at both the real and 
the virtual graphs together. 

A. 5. 2. The transverse part of the ladder diagram 

The calculation of the ladder graph has been carried out as outlined in the main 
text, yielding: 

f- 2 C F du f 1 dz—^—f [1 - 9{z - z cut )} . (A.36) 

8n 2 Jo Jo z(l — u) VI — uj 

The second term, with its ^-function, is the collinear subtraction, whose UV di- 
vergence is cancelled by subtraction of the asymptote. The above formula assumes 
that z cut < 1. If <2 cu t > 1, then we must extend the z integration in the second term 
beyond the limit z — 1, of course. 
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A. 5. 3. Longitudinal part of the vertex correction, self-energy and ladder 
diagram 

The calculation has been carried out as outlined in the main text and yields the 
following for the ladder graph: 



The real and virtual graphs give for both the self-energy and the vertex correction. 
There is no virtual graph for the ladder diagram. 

A. 6. Conclusions 

• We have a systematic method for estimating the sizes of higher-order graphs. 

• The natural expansion parameter in QCD is of the order of 



In answer to a question asked when this work was presented at a conference, let 
us observe: The above number may not be the actual expansion parameter, but 
we argue that this natural size sets a measure of whether actual higher-order 
corrections are of a normal size or are especially large. 

• The method allows an identification of appropriate sizes for renormalization 
and factorization scales. 

• By an examination of the kinematics of graphs, we can identify contributions 
that are large compared with the estimates based merely on the size of the 
natural expansion parameter. 




(A.37) 




x group theory x number of graphs. 



(A.38) 
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• Our method should give a systematic technique of diagnosing the reasons 
for the large corrections, and hence of indicating where one should work out 
resummation methods. 

• The method gives an algorithm for the numerical integration of graphs, both 
real and virtual. 

Of course, there is also a lot of overlap with work such as that of Catani and 
Seymour [ST] , of Brodsky and Lu [53) , and of Neubert p5|. 
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